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I  INTRODU'CT.iON 


The  purport!  of  -lie  work  performed  under  this  contract  was  to 
investigate  experimentally,  the  accuracy  with  which  a  target  moving 
dOo\e  the  ionosphere  can  be  located  an*]  tracked.  The  principal  line 
of  atti  k  on  the  basic  problem  has  been  to  measure.,  by  radio  techniques, 
tne  appa^erj'^  p.)sit.ion  of  radio  satelli+es  as  they  pass  wi-".hin  the  line 
of  -'.c^gnt  o:  a  ground  I’eceiving  station,  an  1  to  compare  the  apparent 
puCLticn  With  the  true  position  over  a  range  of  ele/ation  angles* 

These  meas.^red  deviations  are  then  compared  with  deviations  calculated 
a'.±ng  av  iiiabie  refradive  index  data  which  a'^’e  t' epresentative  of  the 
a/ea  and  time  in  which  the  measurements  were  raken  Since  it  ivS 
ecognized  '.hai  ^'hese  devia'^.ions  are  greatest  at  low  elevation  angles, 
and  '..na'  1"^  i:  1  g:  ea*^  esi.  iinpcr  !.an( .e  ''.o  let.e  "t  and  i  rack  at  low  eleva- 
n  mg^e- ,  "he  w  ' K  under  this  contract  has  been  -onrentrated  on 

i  ,  .  _o 

i.'  n  ingj^e:  Dexcw  . 

The  bi'c  approa-^h  been  '  o  use  ine  cea-^n'^  er feromet er  technique 
:  mca^'.. .  u.g  hau  epp^r'ert  elevation  angle  of  the  saieilite  versus  time 

Ci  '^ng  tti-’n  Pa'"'’.  hi  -'rcer  t.j  be  ab'U  to  determine  tne  true  trajectory 
'  riteVjd-  r  u)gie  tj.me),  during  a  gi'ven  pass,  simultaneous  doppJer 

•  rr.-i  ij.c.'e  ji  i  b^.cn  ^o'clnei.,  The  te’Ctnique".  for  using  doppler  data  for 
'  zC.\  Ci.rjg  ’he  r.^L-  i  je.  t'V,  luc  i:  OirA/assed  in  Gome  detail  in  Section 
V  By  mp  1  :  ng,  ippj"''en’  -in  i  'rue  eieva'^icr  angle  as  the  :atellite  rises 
’'ne  d)*  iin  a  u:'’ve  jT  ^lie  elev-r.ion  angle  er'oi  versus 
'he  appiver*^  ingie, 

Oro-  ’f  Uic  j].;i  -  f  ^nis  ;:orb-»’  3'  '  is  *•  c.  de^  eTmiine  now  well  one  an 
•  mn*.  i  .ACt  '  i  r  A  a  f-’cm  avrij  i  iblo  icfra  •tive  index  data. 


AcoOiTviingiy ,  the  measured  and  computed  ^  curves  are  compared  under 
various  refraction  conditions,  and  conclusions  are  drawn  concerning  the 
uncertainties  in  making  an^lar  refraction  correotionSe 

It  also  ha?  been  found  possible  to  relate  radar  range  uncertainties 
tc  the  observed  fluctuations  in  doppler  frequencies,  thus  enabling  one 
to  estimate  the  degree  of  reliability  with  which  range  error  corrections 
can  be  made. 

It  was  found  during  the  investigations,  that  signals  could  be 
received  from  the  satellite  while  it  was  still  several  hundred  miles 
below  the  norizon*  These  "precursor”  signals  represent  one  of  the  un¬ 
expected  bonuses  of  the  work  done  under  this  contract <>  They  indicate 
that  one  can  expect  to  detect  missiles,  for  example,  well  before  they 
come  up  over  ^  he  horizon.  Because  of  the  importance  of  this  phenomenon, 
^ne  c harac e’^i'^tics  of  precursors  have  been  analyzed  in  more  detail, 
aitnougn  hic  ^ork  was  not  originally  called  fo^’^.  Recommendations  have 
been  made  for  further  investigating  the  precursor  mechanism  in  order  to 
exploit  i-^.?  use  for  early  letection 

Final xy,  one  can  estimate  the  maximum  degree  of  uncertainty  in 
targe^  jo.'i‘'iwn  ausci  by  the  very  rapid  fluctuations  in  Such  flue- 
tua*^iori.  1  the  depth  of  observed  interference  minima  because  they 
''e^ui+  in  a  targe'*  having  an  apparently  finite  sise<,  Vfhile  other  pro- 
i'esrer  tend  to  fill  in  the  minima,  one  can  determine  apparent 

t  i*^get  dimjenoions ;  fa.  which  :an  be  considered  as  a  safe  upper  limit, 
on  ■*ne  Dari^  oi  maxim-om  observed  depth  of  '^.he  minima. 

In  tne  pro"ev^'-'  of  performing  xhe  experimental  studies,  it  has  been 
necessary  **0  bi^lld  a  versatile  field  facility,  including  large  tracking 
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antennas,  narrow  band  phase  lock  receivers,  (capable  of  tracking  ioppler 
frequency  changes) ,  internal  timing  standards  which  need  to  be  locked 
to  WWV  only  once  a  day,  and  multi-channel  recording  equipment  for 
simultaneously  record:^ ng  interferometer  data  on  several  frequencies. 

The  equipment  details  are  described  in  Section  III,  and  the  general 
description  of  the  Point  Buchon  Field  Site  is  given  in  Section  11. 

As  a  result  of  the  work  performed  to  date,  much  valuable  infor¬ 
mation  has  been  obtained,  and  avenues  have  been  indicated  for  obtain¬ 
ing  much  new  information  for  providing  answers  to  basic  questions  on 
target  location  accuracy  and  early  detection «  The  equipment  complex 
for  performing  the  past  work  can  be  readily  expanded  to  perform  the 
new  work  which  is  recommended . 
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II  PQIOT  BUCHQN  FIELD  SITE 


The  experimeiivaJ  measurements  described  in  this  report  were  carried 
out  at  the  3BA  field  site  located  on  the  California  coast  near  Point 
Buchonc  Figure  1  shows  the  location,  about  170  miles  northwest  of 
Los  AngeieSo  In  selecting  this  Site,  several  requirements  were  satis- 
fieOo  First.,  it  has  a  wide  (170°)  miimpeded  horizono  This  is  essential 
for  sea-inter  fee ometer  measurements  on  a  maximum  number  of  satellite 
passes,  Of  pa^^ticula'^  advantage  at  this  site  is  a  clear  horizon  some¬ 
what  east  of  due  south  (to  160°  azimuth)  which  allows  measurement  of 
a  n'omber  of  satellite  sets  as  well  as  rises «  A  recond  requirement 
for  a  sea-interferometer  site  is  a  steep  foreground  slope.  The  Point 
Buchon  site  has  a  maximum  depression  angle  of  19°  to  the  water^s  edge* 

The  depression  angle  exceeds  10°  over  the  range  from  173°  -  285°  azi¬ 
muth  «  Figu-’e  2  chows  ^he  maximum  depression  angle  in  mllliradians 
versus  azimuth  for  the  Highest  antenna  site^  Since  we  are  concerned 
with  meacuring  at  eievation  angles  below  10°  ,  tnece  allowaole  depres¬ 
sion  angJei  have  been  adequate  for  observing  most  visible  passes 
occurring  in  ^ne  we:  tu^^n  naif  of  the  sky* 

A  ’  nira  '.ar  t  requirement  which  is  offered  is  a  low  radio  noise 

environment c  The  site  is  located  on  a  private  ranch  and  has  very  little 
automobile  ’^*’afti  *  Other  sources  of  noise  are  hidden  by  mountains  to 
she  ea^t.  The  cniy  interference  problem  has  been  occasional  radar  noise 
from  the  Vandei  bu^-g  Air  Force  BavSe  area,  50  miles  to  the  southeast*  This 
has  interierei  WLtn  only  three  out  of  nearly  70  passe*^  recorded  to  date* 
The  a/ailabiiity  of  commercial  power  nearby  is  also  a  distinct  asset, 
since  affora.r  mum  better  frequency  and  voltage  regulation  than  is 
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possible  uith  a  portable  power  plant „ 

Figure  3  shows  a  profile  o^  the  slope  on  which  four  antennas  are 
installed^  at  locations  marked  A,  B,  C,  and  D.  The  elevations  of  the 
antennas  are  approximately  270  feet,  30C  feet,  760  feet  and  860  feet  msl, 
respectively,  Tne  lower  two  antennas  (A  and  B)  are  9'  x  parabolic 
sectors  separated  by  about  150^  horizontally*  An  instrrmient  trailer  is 
located  near  the  upper  of  the  two*  This  area  is  just  below  an  existing 
ranch  road  and  the  commercial  powerline*  The  upper  pair  of  antennas 
(C  and  D)  are  18'  x  45'  parabolic  sectors  and  are  separated  by  about 
250'  horizontally*  A  second  instrument  trailer  is  located  near  the 
upper  antenna  of  this  pair  also.  The  upper  area  is  about  500  feet 
higher  and  1400  feet  northeast  of  the  lower  area*  Power  is  delivered 
to  this  area  via  a  ±500  foot  run  of  parkway  cable  ^3-^2  conductors) 
laid  on  tile  ground  to  reduce  poverline  noise*  Power  is  supplied  to 
both  areas  at  480  ^  -  3  phase  for  operating  the  antenna  drive  motors* 

A  portion  is  transformed  do^wn  to  120v  single  phase  for  the  electronics 
equipment  and  ether  needs* 

Both  instrument  trailers  are  18  feet  x  8  feet  insulated  aluminum 
shells  whi  h  were  made  by  Masterbuilt  Trailer  Company  in  Los  Angeles. 

They  were  ‘^.-^w’ed  to  tnei?"  respective  locations,  the  wheels  were  removed 
and  ea^'h  '^nell  was  levelled  and  supported  on  niomerous  wooden  blocks* 

Each  unit  has  a  thermostatically  controlled  fan  for  exhausting  air 
heated  by  dissipation  from  the  electronic  equipment*  When  the  inside 
air  temperature  ^  icec  above  a  set  level  ,  a  vane  is  automatically  moved 
to  direct  ■'.ne  iieated  air  outside*  In  this  way  the  electronics  equip¬ 
ment  IS  kept  at  a  reasonably  constant  temperature,  which  is  important  for 
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maintaining  '.he  frequen::'y  stability  of  the  various  receivers,  coiinters, 
etc. 

Figi^re  4  shows  the  C-D  area.  The  two  antennas  are  shown,  together 
'With  the  instiaimen'-'  trailer  at  the  D  citco  A  portion  of  the  half- 
mile  road  which  was  constructed  to  the  area  is  also  shown.  On  the 
elope  just  above  the  trailer  may  be  seen  a  radiosonde  receiver, 
which  is  metalled  for  obtaining  on-site  refractive  index  data  for 
the  troposphere. 

The  purpose  in  laying  out  the  antennas  in  pairs  is  to  provide  two 
separate  sea-interferometers  at  each  of  two  frequency  ranges.  Separat¬ 
ing  the  heights  of  each  antenna  by  about  10^,  results  in  the  interference 
minima  in  one  pattern  being  displaced  significantly  from  those  in  the 
companion  pattern.  This  serves  the  double  purpose  of  helping  to  elimi¬ 
nate  pattern  ambigui'^ies  as  well  as  providing  twice  as  many  points  where 
apparent  elevation  is  accurately  determined.  The  patterns  of  the  upper 
pair  of  an'^enna:  h^ve  minimn  spacings  of  about  6mr  at  108mc  and  3nir 
ac  2l6mc.  Experience  to  date  has  shown  that  the  upper  frequency  limit 
for  these  anienna*"  ic  probably  about  400mc.,  provided  simultaneous  lower 
frequency  data  are  taken  to  resolve  pattern  ambiguities. 

The  luv-^er  pair  of  antennas,  at  roughly  one-third  the  height  of 
t.he  appear  pai"^.  have  interferometer  patterns  three  times  as  broad,  and 
■therefore  can  be  used  up  to  frequencies  of  the  order  of  12C0mc„ 

Measuiements  to  date  have  been  in  the  ^>4-324  me  range,  such  that 
the  upper  antenna  pair  has  been  c-idequate..  Future  work  on  higher  fre¬ 
quencies^  Will  ne'.ewcsitate  operation  of  the  lower  antenna  pair.  When 
measurements  cover  a  range  of  frequencies  above  and  below  400  me,  both 
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paiTo  will  be  operated  simultaneously , 


TIT.  DESCRIPTION  OF  EQUIPMENT 

The  field  site  instrumentation  at  Point  Buchon  consists  of  an  integrated 
multireceiver  installation  vrith  associated  antennas,  recorders,  timing  cir¬ 
cuits  and  frequency  standards  for  measuring  and  recording  the  apparent 
elevation  angle  of  satellite  radio  sources  by  sea-interferometer  techni- 
•pies.  The  overall  block  diagram  is  indicated  below* 
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Each  f  the  major  components  of  the  instrumentation  is  described 
below. 


Antenna:' 

Figi'^e  :  1C  a  closeup  view  of  the  D  antenna.  Details  of  the  mounting, 

-oU  h  ar  +he  '  on^.rete  block  counterweights  are  illustrated,  together  with 
the  25-foot  long  feed  boom  assembly.  This  assembly  is  entirely  non-metallic . 
Each  boom  is  a  hollow  tapered  structure  made  from  glued  plywcod  and  covered 
wi’^.h  a  1/8"  sheatn  of  fiber  glass  and  resin.  The  assembly  is  supported 
by  fiber  glass  ropec  All  metallic  fittings  are  located  either  at  the  re- 
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fleeter  face  or  at  the  feed  structure  end  of  the  assembly.  The  feed 
structure  has  five  horizontally  polarized  arrays  and  one  vertically 
polarized  feed.  Both  the  C  and  B  antennas  can  receive  simultaneously 
on  54,  108,  162,  216,  and  324  me. 

The  reflectors  are  constructed  of  aluminum  tubing  and  faced  with 
expanded  aluminum  mesh.  Each  structure  was  designed  to  withstand 
40mph  winds  unaided.  During  nonoperating  periods,  each  antenna  is 
secured  to  a  pair  of  vertical  tubular  steel  masts,  which  in  turn  are 
guyed  by  steel  cables.  Vhen  secured  in  this  way,  the  antennas  can 
withstand  winds  in  excess  of  100  mph. 

Both  the  C  and  D  antennas  are  fully  steerable  in  elevation  and 
azimuth.  The  mounting  was  designed  to  allow  continuous  180^  of 
elevation  tracking.  In  this  way,  zenith  passes  may  be  tracked  without 
the  neceSwSity  of  rotating  180^  in  azimuth  when  the  satellite  passes 
through  the  zenith. 

The  1  h.p,  drive  motors  are  capable  of  turning  speeds  of  2  rpm 
in  azimuth  and  1  rpm  in  elevation,  although  satellite  tracking  speeds 
average  much  slower.  The  half-power  beam  widths  are  about  14°  (hori¬ 
zontal)  by  35°  (vertical)  at  108mc  and  about  half  of  this  at  2l6mc . 
Satellite  motion  in  the  lowest  10°  of  elevation  is  such  that  one  can 
keep  the  satellite  in  the  beam  by  shifting  azimuth  at  about  one  minute 
or  half-minute  intervals.  To  obtain  doppler  data  over  the  entire  pass 
requires  more  rapid  motion,  particularly  at  nearest  approach,  but  step¬ 
wise  tracking  is  still  quite  satisfactory  in  this  frequency  range. 

The  antenna  controls  were  designed  to  supply  power  to  the  drive 
motors  at  eijther  120  v,  240  v,  or  4B0  v,  depending  upon  the  torque 
requirements .  Sels^m  indicators  are  provided  in  the  equipment  trailer 
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for  remotely  indicating  azimuth  and  elevation  angles  of  both  antennas. 
In  practice,  an  antenna  turning  schedule  is  prepared  in  advance  of 
each  pass,  indicating  the  azimuth  and  elevation  angles  to  be  maintained 
at  each  time  when  a  change  is  required.  The  operator  then  uses  this 
schedule  for  both  the  C  and  D  antennas;  changing  to  the  new  positions 
at  the  times  indicated.  Since  the  operator's  post  is  adjacent  to  the 
data  recorder,  verbal  instructions  can  readily  be  given  if  it  appears 
that  the  satellite  is  drifting  out  of  the  antenna  beam. 

Receivers 

Twelve  receivers  make  up  the  receiver  complement  of  the  Point 
Buchon  I'ield  Site^  These  receivers  are  SRA  Modular  Receiving  Systems. 
In  the  configuration  used  fixed  tuned  rf  amplifier  modules  and  crystal 
v.ontrolled  local  oscillator  modules  are  used  in  the  receiving  systems 
to  allow  maxim^im  sensitivity  to  be  realized  for  the  reception  of  the 
signalo  from  each  satellite  tracked.  The  receivers  are  used  in  either 
^f  two  modes  of  operation:  (l)  Fixed  tuned,  or  (2)  phase  lock  fre¬ 
quency  tracking.  The  design  and  construction  of  the  SRA  Modular 
Receiving  Systems  allow  the  required  llinctionally  separate  units  of 
the  receiver  to  be  plugged  in  to  the  receiver  shell  so  that  all  the 
characteristics  of  the  receiving  system  can  be  chosen  by  appropriate 
selection  of  the  particular  type  units  to  be  used.  A  block  and  module 
diagram  of  one  of  the  phase  locked  frequency  tracking  receiving  systems 
IS  shown  below,.  This  receiver  has  the  following  characteristics: 


1, 

Frequency  of  operation 

]07.940  to  108.090  me. 

2, 

First  IF  Frequency 

20.500  mco 

Second  IF  Frequency 

455  kcs. 

4. 

Pre-detection  bandwidth 

500  cps. 
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Noise  figure 


2o5  db. 


6o  Phase  locked  tuning  range  +  5  kc « 

7n  Phase  error  at  20  db  S/N  ratio;  15^  maximum  at  50 


Measurements  of  the  frequency  of  the  received  signal  as  well  as  the 
amplitude  of  the  "'eceived  signal  are  recorded  at  the  outputs  of  this 
receiver c  In  peration,  the  receiver  tuning  tracks  the  instantaneous 
frequency  of  the  received  signal «  Receivers  ^i^uth  similar  characteristics 
it  operating  frequencies  around  l62  mcs,  2l6  mcs,  and  3^4  indc  are  assembled 
when  different  satellites  are  to  be  measured,  merely  by  replacing  one  or 
two  modules  in  the  receiving  system «  The  modular  design  concept  has 
made  p^a^"'  cal  the  taking  of  measurements  on  two  satellites  of  different 
operating  frequencies  such  as  Tiros  and  Transit  with  as  little  as  10 
minutes  available  for  equipment  change-over c  The  operating  frequency 
of  five  receivers  are  changed  without  sacrifice  of  the  optimum  sensiti¬ 
vity  or  selectivity  available  in  receivers  with  performance  peaked  at 
the  opera‘s  ing  frequency^ 
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The  fixed  tuned  crystal  controlled  receivers  are  assembled  vdth 
many  of  the  ^ame  modules  a?  the  phase  lock  frequency  tracking  receivers 
described  above «  The  455  kcs  second  amplifier  employs  a  6  kc  filter 
which  aJlows  the  doppler  shift  expected  in  the  received  satellite  signals 
to  be  contained  in  the  wider  bandwidth  system «  A  block  and  module  diagram 
of  this  t,ype  receiver  is  the  same  as  that  shown  above  with  the  exception 
that  the  reference  OvScillator,  phase  comparator,  and  loop  filter  are  not 
employed* 

Frequency  Measuring  System 

The  frequency  of  the  received  signal  is  measured  by  electronically 
counting  the  operating  frequencies  of  the  local  oscillator  voltages 
in  the  super^hetrodyne  receiver*  Recordings  of  the  doppler  frequency 
made  by  digital  printing  of  the  counted  frequency  of  the  1st  local 
oscillator  wnich  in  turn  is  phase  locked  to  the  received  signal.  The 
electronic  counter  used  is  the  Hewlett  Packard  5240.,  shown  at  lower  left, 
Figure  6c 

Recorders 

Two  recording  systems  are  used  at  Point  Buchon  to  collect  the 
satellite  data;  ^j)  an  eight  channel  ink-pen  recorder  (Offner  type 
RC  bynograph) ,  shown  at  right  center  of  Figure  6.  Each  of  detected 
recei'ver  outpuls,  which  arej  proportional  to  input  signal  strength  are 
simultaneously  recorded  as  analog  traces  on  a  16”  recorder  paper; 

(2)  a  digital  p-^inter  which  prints  the  frequency  of  a  local  oscillator 
pnase  locked  to  one  of  the  oignals  received  from  the  satellite* 

Gide  pens  are  used  on  the  pen  recorder  to  make  second  and  minute 
intervals  on  ^he  paper  chart.  The  time  is  printed  directly  on  the 
'ligital  doppJer  frequency  record.  Samples  of  thewse  records  are  shown 


in  Figure 


Timing  Circuitn 

A  cryctal  controlled  secondary  frequency  standard  is  incorporated 
VI thin  the  electronic  counter,  and  normally  supplies  the  gating  commands 
to  the  electronic  timer o  An  output  is  taken  from  this  frequency  standard 
and  divided  until  accurate  one  second  and  one  minute  Impulses  are  avail¬ 
able  „  Controls  aie  available  to  vary  the  frequency  divider  ratios  in 
order  that  the  timing  pulses  at  the  output  of  the  time  scales  can  be 
brougnt  into  synchronism  with  WWV  standard  transmissions «  These  timing 
circuits  control  the  timing  of  the  electronic  counter  gate,  the  second 
and  minute  marks  on  the  recorded  data,  and  synchronizing  signals  to  the 
shutter  on  the  ballistic  camera o  The  equipment  is  shown  at  upper  left, 
Fig-ure  6o 

The  Geo Dhyc leal  Sensors 

ho^;eT•ol  instr-oments  are  available  at  Point  Buchon  to  measure  and 
mcni'^.Gr  ^ne  geophyric-al  environment  at  the  site*  Among  these  are: 

i*  Tide  gauge,  2*  Ballistic  camera,  3-  Radiosonde* 


Tide  Gau^.e 

The  "ide  gauge  is  in  the  form  of  an  air  chamber  with  a  filtered 
jpening  at  tne  bottom  and  a  pressure  connection  at  the  top*  The  unit 
i:  bu-riea  seve""al  feet  below  the  sand  level  at  a  point  on  the  beach 
above  mean  tide  level  *  The  beach  sand  serves  to  filter  out  individual 
Wives,  but  allow:  the  averaged  sea  level  at  a  given  time  to  be  estab¬ 
lished  above  tne  cnamber*  The  pressure  inside  the  chamber  is  trans¬ 
mitted  via  j /a”  copper  tubing  to  a  pressure  transducer  located  well 
above  high  tide.  This  unit  is  a  potentiometer  whose  sliding  contact 
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IS  controlled  by  the  input  pressure*  It  has  a  total  range  of  C-5  psi, 
corresponding  to  a  10  foot  tidal  range;,  which  is  adequate  for  the  Point 
Buchon  coastline*  The  transducer  is  part  of  a  battery  operated  bridge 
circuit,  wnose  unbalance  voltage  is  nearly  linear  with  the  head  of 
water  above  the  air  chamber  inlet.  The  output  is  recorded  on  a  0-1 
raa  Esterline-Ang^is  recorder,  with  a  full-scale  sensitivity  of  0-10 
feet  of  sea  water.  The  bridge,  battery  and  recorder  are  housed  in  a 
wdtertighr.  box  on  the  edge  of  the  cliff  overlocking  the  beach  where 
the  gauge  is  located.  During  operations,  the  chart  is  checked  daily 
and  picked  up  at  tne  end  of  each  operating  period. 

An  attempt  was  made  to  install  the  air  chamber  by  burying  it  a 
foot  below  the  beach  level  at  the  water’s  edge  at  low  tide.  This 
depth  of  burial  was  insufficient,  and  wave  action  uncovered  the 
chamber  a  few  nours  af^er  it  was  installed.  During  the  operating 
in*errii,  noweve''^ ,  it  operated  as  expected,  showing  that  the  principle 
1.-  a  fearibxfc  one.  A  new  chamber  has  been  constructed  which  can  be 
driven  into  the  sand  several  feet.  It  is  planned  to  install  it  fur¬ 
ther  bark  on  the  beach  at  a  point  where  5  feet  burial  will  place  it 
just  below  ne  lowest  tide  level.  No  further  difficulties  are  expected 
in  dide  mearuremen’:.. 

Ballisbi'.  Camera 

This  IS  a  fixed-focus  camera  with  a  20”  e.f.l.F-5r6  serial  camera 
len:  „  It  ha->  plate  Holders  for  8”  x  10”  or  4”  x  5”  cut  film.  Using 
the  linger  fi im  gives  a  field  of  view  of  20°  x  25^.  The  camera  has  a 
ruuating  paddle  ju.  t  in  front  of  the  lens,  which  is  driven  by  a  1  rps 
'ynThronoUw  motor  whj.ch  is  synchronised  with  the  timing  ticks  pro¬ 
vided  as  described  above.  The  paddle  nerves  to  interrupt  the  satellite 
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track  at  known  one-second  intervals ♦  The  complete  unit  is  mounted  on  a 
rugged  tripod,  which  al  lov/s  it  to  be  clamped  in  any  desired  position, 
for  photographing  various  portions  of  satellite  trajectories.  Using 
Royal  Pan  film,  the  camera  records  stars  down  to  7th  magnitude.  With 
fast  moving  satellites  however,  5th  magnitude  (barely  visible)  will 
be  required  to  obtain  usable  tracks. 

Radiosonde 

The  radiosonde  equipment  installed  at  Point  Buchon  is  the  AN/GMD-IA 
rawinsonde  wldch  automatically  tracks  a  balloon-borne  transmitter^  The 
ground  equipment  can  be  seen  on  the  slope  above  the  D  antenna  location  in 
Figure  4o  A  data  recorder  is  located  inside  the  instrument  trailer «  It 
records  temperature,  pressure  and  humidity  data  transmitted  by  the  balloon- 
borne  transmit ler„  These  \d.ll  be  launched  at  low  altitude  northwest  of 
the  receiver  location,  such  that  they  can  be  tracked  during  all  parts 
of  their  ascent. 
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TV  .  EXPERIMENTAL  DATA 


The  basic  radio  data  which  are  obtained  during  a  given  satellite 
pass  are  in  two  for.ns:  (l)  a  multi-channel  graphic  recorder  chart  and 
(2)  a  digital  print-out  paper  t^pe.  The  interferometer  data  consist 
of  simultaneous  traces  of  received  signal  level  versus  time  for  each 
receiver  e  Nomally,  two  receivers  are  used  for  each  frequency;  one 
receiver  using  the  C  antenna  and  the  other  using  the  D  antenna «  On 
Tiros  passes  only  one  frequency  (l08mc)  was  transmitted,  so  only  two 
interferoinett  r  traces  are  obtained  on  these  passeso  Most  Transit 
passes  have  at  least  four  traces,  two  on  2l6  me  and  two  on  l62  mCc 
The  54-324  me  transmitters  were  not  operating  during  most  passes,  but 
one  trace  each  was  obtained  on  these  frequencies  as  well,  when  possible. 

As  described  in  the  previous  section,  one  receiver  on  each  pass 
was  a  phase-loek  type,  with  500  cps  effective  bandwidth.  Its  signal 
output  was  recorded  as  one  of  the  interferometer  traces,  and  in  addi¬ 
tion,  the  local  oscillator  frequency  was  recorded  as  a  measure  of  the 
doppler  frequency*.  One  channel  of  the  graphic  recorder  was  used  to 
record  the  doppler  graphically,  and  in  addition,  the  L.O.  frequency 
was  printed  out  digitally,  along  with  time,  at  two  second  intervals  on 
a  -eparate  paper  tape^ 

Timing  infoimation  for  the  digital  tape  as  well  as  for  time  ticks 
on  the  multi-channel  graphic  chart,  were  derived  from  a  common  timing 
source  to  preserve  accurate  time  correlation  between  the  two  forms  of 
recorded  data.  Time  ticks  on  the  graphic  chart  were  at  1  minute  inter¬ 
vals  along  one  margin  and  at  1  second  intervals  along  the  other  margin. 
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A  chart  speed  of  5  nun  per  second  was  used  during  the  pass,  allowing 
time  interpolation  to  +  0.1  sec.  Figure  7  shows  a  recording  of  a  ^ 

typical  Transit  pass*  The  upper  part  of  the  figure  shows  the  multi¬ 
channel  graphic  chart  and  the  lower  portion  shows  a  section  of  the 
digital  aoppler  data« 

In  addition  to  the  radio  data  described  above,  one  needs  other 
data  j.n  order  to  work  up  each  pasSo 

The  exact  height  of  each  antenna  above  sea  level  is  required  in 
order  to  assign  apparent  elevation  angles  to  each  of  the  interference 
minima o  The  height  of  each  antenna  above  mean  sea  level  was  determined 
by  vSurveys  conducted  by  a  team  from  RADC.  This  team  installed  three 
bronze  survey  markers ,  at  locations  near  each  antehnasite.  The  geo¬ 
graphic  co-ordinates  as  well  as  the  height  above  a  local  bench  mark 
were  determined  by  the  team  for  each  of  these  markers «  After  installa¬ 
tion  of  the  antennas,  these  markers  were  used  to  determine  the  co¬ 
ordinates  and  heights  of  the  intersection  of  rotation  axes  for  each 
antenna « 

The  heights  above  mean  sea  level,  as  determined  above,  are  then 
corrected  for  tlie  tide  level  during  each  pasSo  Tide  tables  have  been 
used  for  this  pu^^pose,  and  have  been  found  to  be  adequate  for  the  +  0,3 
foot  accuracy  required c  The  actual  tide  level  was  checked  by  a  survey 
to  tne  Water*-  edge  and  was  found  in  agreement  with  the  tide  table  to 
within  Ool  feeto  It  is  realized  that  occasional  deviations  will  occur 
due  to  offshore  storms,  and  an  attempt  was  made  to  install  a  local  tide 
gauge.  Although  the  attempt  was  unsuccessful  as  described  in  Section  III, 
it  ivS  planned  to  complete  the  installation  in  the  near  future.  Meanwhile, 
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the  tide  table  corrections  have  been  found  to  be  sufficiently  reliable 
for  past  purposes o 

One  type  of  auxiliary  data  is  needed  prior  to  each  pass.  This  is 
in  the  form  of  predicted  satellite  trajectories,  which  are  required  for 
prograrmning  the  antenna  pointing  angles.  Since  numerous  passes  occur 
during  any  given  day,  and  since  the  usual  observing  periods  extend  over 
about  two  weeks y  a  method  was  devised  which  allows  the  field  team  to 
pxepare  its  om  pointing  schedule  for  those  passes  on  which  it  is 
desired  to  mL,ke  measurements. 

First  y  a  senes  of  general  charts  was  prepared  for  each  satellite 
to  be  utilized.  Figure  8  is  an  example  of  one  such  chart.  The  family 
cf  curves  with  arrow  indicates  a  series  of  typical  tra jectories(eleva- 
tion  angle  vs  azimuth)  which  the  satellite  wil]  follow.  A  circular 
orbit  [lero  eccentricity)  was  assumed  for  these  curves,  since  the  eccen¬ 
tricities  of  tne  satellites  of  interest  are  relatively  small.  A  second 
set  of  Curves  is  superimposed  on  the  trajectory  curves.  The  intersections 
of  these  curves  with  a  given  trajectory  curve  indicates  the  expected 
position  of  tne  sa'^.ellite  at  one  minute  intervals. 

In  order  to  specify  the  trajectory  of  the  satellite  during  a  given 

pass  using  tnese  shirts,  one  only  needs  to  know  the  predicxed  time, 

elevation  angle  and  azimuth  of  one  point  during  the  pass.  This  point 
IS  plotted  on  the  appi’opriate  chart,  and  the  remainder  of  the  trajectory 
can  be  easily  sketched  in  by  interpolation  between  adjacent  curves  on 

the  chact.  Knowing  the  time  of  arrival  of  the  satellite  at  the  given 

point  allows  one  to  estimate  the  times  of  arrival  at  all  other  points 
along  the  pasr.  If  the  satellite  rises  earlier  than  predicted,  the 
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estimated  trajectory  is  moved  eastward  by  an  appropriate  amount, 
knowing  that  the  time  intervals  between  adjacent  trajectory  curves  is 
20  minutes.  Similarly,  if  the  satellite  rises  later  than  expected, 
the  estimated  trajectory  is  moved  westward,  to  account  for  the  earth* s 
rotation  during  the  interval. 

In  order  to  preserve  the  charts  for  repeated  use,  it  has  been 
found,  desirable  to  use  an  overlay  of  transparent  plastic  film,  such  as 
Saran  Wrap,  and  plotting  the  trajectory  with  a  grease  pencil.  Since 
Saran  Wrap  does  not  adhere  well  to  paper,  it  n s  recommended  that  the 
charts  themselves  be  laminated  in  plastic.  In  this  way,  the  overlay 
film  will  adhere  to  the  chart  without  the  necessity  of  taping. 

Witxi  the  aid  of  a  set  of  these  charts  (one  for  northbound  passes 
and  one  for  southbound  passes)  for  each  satellite,  the  pass  predictions 
need  only  include  one  point  (usually  near  rise)  for  each  visible  pass 
occurring  during  the  two  week  period  of  interest r  These  predictions  are 
computed  using  Space  Track  orbital  data  and  an  electronic  computer. 

Past  experience  with  Tiros  and  Transit  satellites  has  indicated 
that  the  method  described  above  is  adequate  for  antenna  pointing  pur¬ 
poses.  During  the  pass,  it  is  occasionally  apparent  that  the  satellite 
13  getting  a  little  ahead  cr  behind  the  prediction,  but  these  cases 
are  easily  remedied  by  leading  or  lagging  the  schedule  by  a  small  amount. 
By  and  large  >  the  m.ethod  has  been  quite  practical  and  the  accuracy  has 
been  all  that  is  required  for  keeping  the  satellite  in  the  antenna  beam. 
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V.  ^ATA  REDUCTION  TECHNIQUES 


Valid  Sea-Interfarometer  Interval 

The  firr.t  otep  in  the  reduction  of  the  interferometer  traces,  is 
to  determine  the  limits  of  the  usable  datao  The  criterion  in  this 
case  is  the  Tnaxiinum  elevation  angle  (Figure  2)  allowed  by  the  fore¬ 
ground  terrain  at  the  azimuth  of  the  particular  pasSr  For  this  purpose, 
the  antenna  pointing  schedule  is  helpful,  since  it  indicates  azimuth 

timoo  By  also  counting  minima.,  one  can  roughly  determine  the  approxi¬ 
mate  elevation  angle  at  various  times «  When  the  approximate  elevation 
ingle  approaches  within  a  degree  of  the  limit  set  by  Figure  2,  it  is 
lonsiiered  that  the  usable  limit  is  reached ^  Beyond  this  point,  the 
pha-e  of  ihe  reflected  wave  v-all  probably  be  displaced  by  terrain 
reflection c  su^h  that  the  elevation  angles  computed  for  sea-reflection 
will  be  invalid.. 

Timers  ot  O'^.currence  of  Minima 

One  next  determines  the  times  of  occurrence  of  each  interference 
minimum  observed  in  the  valid  interval®  Most  of  the  minima  are  suffi- 
''ion'': 'iy  ^ha'^p  {'’ee  Figure  7)  that  one  can  read  the  times  to  jp.l  second 
(O  5  mm  on  the  -’ecorder  ^hart)  ®  The  minima  which  are  somewhat  indistinct 
recend)  arc  tabulated  with  a  question  mark,  and  those  which  are 
."till  ii  i.c.uvt  0c5  second  or  more)  are  di^^carded.  The  time 

os‘^Mibj.:^shed  ci:  tne  recorder  chart  by  the  one  minute  ticks  on  the 
upper  margin  and  the  ont.  second  ticks  on  the  lower  margin,  is  adequate 
to  establish  tlie  times  oi  occurrence  of  minima  within  the  prescribed 
Ool  se'^ond  tolerance^  This  tolerance  was  set  by  the  average  angular 
'•rite  of  rise,  whi'''h  is  ibout  one  mr®/3econd  and  the  required  angular 
ac^uri'-y  of  0.1  mr 
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Elevation  Angles  of  Minima 


In  order  lu  plot  the  apparent  trajectory  (apparent  elevation  angle  vs 
time)  for  a  particular  pass,  one  needs  to  determine  the  elevation  angle 
corresponding  to  each  of  the  well-defined  minima. 

In  the  sea  interferometer,  the  elevation  angle  corresponding  to  each 
minimum  is  derived  as  follows: 

F^^r  a  transmitter  at  infinite  distance,  the  path  difference  between  a 

direct  wave  and  one  which  is  reflected  from  the  sea  surface  is 

A  ]  2h  sin  a  (l) 

o  o 

where  h^  is  the  height  of  the  receiver  above  a  plane  tangent 
to  the  earth  at  the  reflection  point,  and  is  the  angle  of 
incidence  of  the  reflected  wave  on  the  sea  surface. 

In  general,  the  reflected  wave  undergoes  a  phase  change  cp.  In  order 
for  tne  direct  and  reflected  waves  to  arrive  at  the  antenna  just  out  of 
phase  tnen  we  require  that: 

A  ]  (n  -  ]/2  +  cp/2n)  X  (2) 

where  X  is  the  wavelength  and  n  is  the  order  number  (],2,3,etc.) 
■f  the  minimum. 

In  the  case  of  a  norizontally  polarized  wave,  cp  is  equal  to  7i  radians 
witnin  \fery  narrow  limits,  regardless  of  wavelength  or  angle  of  incidence. 
The  condition  for  an  interference  minimum  with  horizontal  polarization  thus 
becomes : 

A  1  "  n  K 
Thus  we  <  an  writer 

5 in  ■-  nX/2h^  (3) 
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Since  varies  as  the  distance  to  the  reflection  point  moves,  we  would 
like  to  express  antenna  height  in  terms  of  the  height  h  above  the  sea 
crarface  itself.  From  geometry  we  can  write  (for  d«a): 

h  •  h  +  Ah  =  h  +  d^/2a 

C'  ; . 

where  d  is  the  distance  to  the  reflection  point  and  a  is  the 
raaJUj  of  curvature  of  the  sea  surface. 

F:rm  geometry  it  is  evident  that; 
n, /d  -  tan  a 

Combining  tne  above  expressions  gives: 

^in  a  tan  a  “  n\/2d 

O  :/ 

cini 

a  a  [-lan  +  (tan"'  2h/a)  ] 

Eliminating  d  gives  finally; 

?  i 

in  a  tan  a  [-  tan  a  +  (tan  a  +  2h/a)  ]  -  n\/2a  (4) 

This  expression  cannot  be  readily  solved  implicitly  for  a^. 

SiTiall  Angles 

■'f  we  confine  our  attention  to  ci^'s  small  enough  that  sin  a  - 
tan  a  ”  a  »  we  can  then  re-write  (4)  as: 


(5) 


a  -  a  ^  n\/2h  -  \^/8ah 

0  0  ' 

This  quartic  is  most  easily  solved  numerically  using  Nevrton's  Method, 


Large  Angles 

2  2h 

When  a  exceeds  about  30  mr,  tan  a  is  considerably  larger  than  — , 
0  0  a 

such  that  the  square  root  teirni  in  (4)  may  be  adequately  represented  as 
a  tnree-term  expansion: 


ban  a  +  — - 

o  a  tan  a 


2  3 

2a  tan-^  a 


Thus  (4)  becomes: 

2 

sin  “  n\/2h(l  -  h/2a  tan  a^)  (6) 

Since  the  second  term  in  the  denominator  is  small  compared  to  the  first, 
one  can  use  an  approximate  value  of  a^,  In  this  term  to  obtain  values  of 
sin  which  are  accurate  to  +  ,01  mr.  The  solutions  of  (5)  and  (6)  yield 
values  of  a^,  which  is  the  angle  of  incidence  of  the  reflected  wave  on  the 
sea  surface.  It  is  more  convenient  to  convert  this  to  a,  the  angle  above 
the  horizontal  at  the  antenna  itself.  The  conversion  at  small  angles  is: 
a  ■■  2a  -  (a  ^  +  2h/a 

At  large  angles  (greater  than  30  mr)  the  correction  is  small,  such  that  a 
two  term  expansion  ' f  a  square  root  term  similar  to  that  in  (4)  may  be 
used.  In  this  case: 


a  -  a  -  h/aa 
0  '  o 

In  practice,  both  (5)  ana  (6)  corrected  to  a  have  been  tabulated  for  h  - 
860  ft.  and  760  ft,  (the  nominal  heights  of  our  sea-interferometer  antennas, 
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and  for  frequencies  of  108  me,  ]62  me,  and  216  me.  To  eorreet  for  tidal 
changes,  tabulations  are  aetually  made  for  one  foot  intervals  of  antenna 
height  eovering  the  six-foot  maximum  tidal  range  expeeted  at  the  field 
site.  Interpolation  ean  readily  be  made  to  +  0.1  feet,  which  brings  the 
accuracy  of  the  elevation  angles  within  0.1  mr  at  angles  up  to  300  mr 
(]7  .2°) . 


Finite  Transmitter  Distance 

The  above  expressions  were  derived  assuming  parallel  direct  and  re¬ 
flected  rays  beyond  the  reflection  point.  If  the  transmitter  (or  radar 
target)  is  at  a  finite  range  R,  then  two  corrections  need  to  be  made. 
First,  the  elevation  angle  of  the  direct  ray  must  be  lowered  by  the  angle 
.motendei  ax  the  cransmitxer  by  the  inter ferometer  baseline.  This  angle 
1:.*  ; 


In  the  fig’iro,  we  also  note  that  the  direct  ray  is  now  slightly  longer 
than  it  w^mlcl  be  if  it  were  parallel  to  the  reflected  ray.  This  added 
path  length  is.; 

m  2h  ^/R 

o 
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The  path  aifrerence  between  direct  and  reflected  rays  is  now  reduced  by 
which  raises  the  elevation  angles  of  the  minima  by  Z^l/2h^  =  Sub¬ 

tracting  this  effect  from  Aa'  leaves  the  final  correction: 

Aa  -b^/R 

The  height  of  a  typical  satellite  used  in  our  studies  is  about  500  miles. 
This  gives  the  following  values  of  Aa: 

a  Aa 

0  -  0 . 14  mr 

1 0  mr  -  0 . 1 6  mr 

100  mr  -  0.20  mr 

300  mr  -  0.30  mr 

Although  these  corrections  are  small ^  they  should  be  applied  for  accurate 
resul ts . 

It  mould  be  pointed  out  that  a,  the  earth’s  radius  which  is  used  in 
the  above  expresrions  should  correspond  to  the  curvature  of  the  sea  surface 
at  the  latitude  f  the  observing  station.  However,  the  maximum  value  of 
tne  terms  involving  a,  is  always  less  than  5  mill iradians .  Since  the 
maximum  variation  of  a,  over  the  earth’s  surface  is  of  the  order  of  0.4% ♦ 
me  could  use  a  value  anywhere  in  this  range  and  be  accurate  to  within 
02  mr  .  It  turns  out  tnat  the  radius  of  curvature  at  latitude  35^N  is 
very  nearly  equal  to  that  of  a  sphere  having  the  same  volume  as  the  earth, 
namely  6371  Km  This  was  the  value  used  in  our  computations,  and  the  re¬ 
sulting  uncertainty  in  the  elevation  angle  from  this  cource  is  certainly 
much  less  than  0.01  mr . 

Tide  C  riections 

In  order  to  retain  the  required  accuracy  of  +  0.1  mr.  in  the  computed 
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elevation  angle,  of  the  minima,  over  the  largest  range  of  angles  in  the 
observational  data,  one  needs  to  establish  tide  corrections  to  +0.3  feet. 
As  mentioned  earlier,  attempts  have  been  made  (unsuccessfully)  to  install  a 
tide  gauge  near  the  receiving  site.  However,  the  purpose  of  this  gauge 
Wd;^  only  to  detect  sea  level  anomalies  caused  by  offshore  storms  and  high 
winds  ,  Since  theae  anomalies  are  rarely  more  than  a  few  tenths  of  a  foot 
in  amplitude,  it  is  believed  that  the  tide  cables  offer  adequate  correction 
for  our  purposes.  During  the  attempt  to  install  the  gauge,  a  survey  was 
made  from  the  j’cference  benen  mark  to  tne  actual  water  level  at  low  tide, 
and  the  measured  level  agreed  with  the  tide  table  value  to  +  0.1  foot, 
the  limit  of  accuracy  of  the  table  and  of  the  survey.  This  agreement 
confirms  the  datum  correction  for  the  area  (difference  between  mean  sea 
level  and  mean  low  water)  and  also  confirms  our  interpolation  between 

two  estahli^heu  secondary  tide  gauges,  one  about  10  miles  southeast  and  one 
about  10  miles  north  of  the  receiving  site  coastline.  The  source  of  tide 
c  'Tjection  information  which  has  been  used  throughout  is:  "Tide  Tables  - 
Wes’  C‘ a:.t  of  Morth  ana  South  America"  published  by  the  U,  S.  Coast  and 
Geodetic  Survey.,  On  the  basis  of  the  above  considerations,  it  is  believed 
that  tne  lido  corrections  made-  for  tne  passes  observed  to  date  have  been 
adequate  for  maintaining  the  required  0.1  mr  accuracy  of  all  measurements. 

Resolution  of  Amoip.uities 

In  orde”  to  plot  a  trajectory  curve  of  apparent  elevation  angle  versus 
time,  one  mu: t  match  the  elevation  angler  as  computed  above  with  each  of 
the  observed  minima.  In  a  classic  pattern,  where  all  minima  are  recorded 
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in  an  undisturbed  whj,  this  is  no  problem  since  the  first  recorded  minimum 
corresponds  to  n  -  1,  the  second  to  n  *  2,  eic  .  However,  in  the  presence 
of  strong  i-efiaction,  iioturbed  conditionSy  or  cccasionally ,  failure  to 
record  ihe  tiTst  few  minima  for  equipment  reasons,  one  is  frequently  faced 
witii  trie  problem  of  esiablisning  correct  order  numbers  without  having  the 
first  few  minima  as  a  guide. 

In  most  ca,>es,  5:imultaeous  interferometer  traces  are  available  from 
botn  the  C  and  D  antenna:  ,  Since  tne  tieights  are  different  for  these  an¬ 
tennas  860^  for  D  and  760  for  C) ,  the  minima  occur  at  different  elevation 
angles.  It  turn-  out  that  the  9tn  minimum  for  the  D  antenna  should  occur 
about  0.2  mr  lowe-^  {0.2  second  earlier)  than  the  8th  minimtam  for  the  C 
antenna.  Using  tnis  criterion,  one  can  establish  these  minima  and  work 
back  to  tne  earliest  observed  minimum  on  botn  traces- 

On  a  few  passes,  two-antenna  aata  were  not  taken,  but  two-frequency 
lata  sn  tne  same  antenna  were  recorded.  In  these  cases  the  ratio  of  order 
numbers  for  coincident  minima  is  proportional  to  the  frequency  ratios,  For 
example,  if  data  were  taken  on  162  me  and  216  me,  minimum  number  4  on  21 6  me 
w:'uld  be  expected  to  be  coincident  with  minimum  number  3  on  162  me..  Actu¬ 
ally,  in  those  ..ase:  ,  sin-e  refraction  on  162  me  is  somewhat  more  than  on 
216  me,  one  usuall}^  find:-  the  tnira  i62  me  minimum  a  second  or  two  earlier 
tnan  tne  fourth  2'^b  me  minimum  Since  this  method  is  less  certain  than 
tne  two-antenna  method,  it  was  used  only  when  necessary. 

There  were  a  lew  cases  where  extremely  strong  refraction  occurred, 
during  which  the  C  antenna  data  indicated  a  few  tenths  of  a  milliradian 
more  refraction  tnan  me  higher  D  antenna.  In  these  cases,  the  ambiguities 
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were  resolved  by  fiisl  establishing  the  true  trajectory  and  then  assigning 
an  order  number  which  gave  nearly  the  expected  refraction  at  the  highest 
e^^evaticn  angles.  Since  the  minima  at  the  highest  frequency  (usually  216  me) 
are  about  2.7  irir  apart,  and  since  one  knows  the  expected  refraction  at  high 
angles  to  well  within  this  tolerance,  resolution  of  pattern  ambiguity 
could  always  be  accomplished  using  this  means  if  no  others  were  available. 

True  Trajectory  Cemputation 

Having  established  the  apparent  trajectory  of  the  satellite  within 
0.1  mr  accuracy  dui'ing  each  pass,  one  must  now  determine  the  true  trajec¬ 
tory  in  order  to  determine  the  elevation  angle  errors. 

Since  the  satellite  is  follovring  an  orbit  around  the  earth,  one  could 
■rnmpute  the  true  elevation  angle  versus  time,  if  he  knew  the  true  orbital 
element.::  whicn  apply  to  each  satellite  pass.  It  is  recognized  that  some 
“f  the  orbital  elements  isruied  by  Space  Track  and  other  agencies  are  not 
sufficiently  precise  to  establish  the  f»"ajectory  within  the  required  0.1  mr 
accuracy.  One  of  the  most  crucial  elements,  and  the  one  subject  to  the 
greatest  secular  variations,  is  the  actual  position  of  the  satellite  in 
it,,  orbit,  If  this  position  departs  by  more  than  0.1  second  from  its 
given  value,  the  "’esulting  trajectory  will  not  be  as  accurate  as  required. 

Thiv":  difficulty  can  be  resolved  if  one  has  a  sufficiently  accurate 
method  nf  fixing  the  time  of  passage  of  the  satellite  through  one  known 
point  during  tiie  pass.  If  the  remaining  orbital  elements  applying  to  the 
pass  were  known  witn  sufficient  accuracy,  one  could  establish  the  time 
element  by  the  measured  fix  point. 
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Doppler  Fix  Method 


Since  it  was  lealized  that  most  passes  would  occur  at  times  unfavor¬ 
able  for  optical  observation  of  the  satellite,  and  that  many  satellites  are 
too  faint  to  observe  with  ballistic  cameras  or  theodolites,  even  under 
'ptimum  lighting  conditions,  a  radio  fix  method  was  felt  to  be  essential. 

The  method  developed  for  the  present  applications  involved  measuring  the 
Joppl er  frequency  shift  of  the  received  radio  signal  over  as  large  a  portion 
of  the  pass  as  possible.  The  time  of  passage  through  the  point  of  nearest 
approach  can  be  established  within  10  seconds,  by  determining  the  time  of 
greatest  time  rate  of  doppler  change.  One  then  updates  the  given  orbital 
elements  to  this  time  and  uses  them  to  compute  a  curve  of  doppler  shift 
versus  time.  The  shift  in  frequency  is  computed  relative  to  the  frequency 
at  nearest  approach,  which  is  still  unknown.  The  time  scale  of  the  compu¬ 
tation  is  referred  to  the  estimated  time  of  nearest  approach,  which  is 
Known  only  to  the  nearest  few  seconds.  However,  if  the  orbital  elements 
used  in  the  computation  are  correct,  the  measured  and  computed  doppler  curves 
will  have  the  same  shape  and  by  suitable  shifting  of  axes,  can  be  made  to 
■oinciae  over  all  but  the  lowest  50  mr  of  the  pass,  where  refraction  dis¬ 
places  the  measured  curve  by  a  few  cps.  If  the  two  curves  do  coincide 
over  the  remainder  of  tne  pass  (usually  at  least  10  minutes  long),  one  can 
be  asv^ured  that  ihe  orbital  elements  are  representative  of  the  pass,  and 
t-hat  tne  true  elevation  angles  calculated  using  these  elements  will  be 
va] id  for  the  pass 

Matching  the  twc  curves  in  thi:>  way,  allows  one  to  determine  the 
time  of  passage  of  the  satellite  through  the  point  of  nearest  approach  to 
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+  0.1  second  find  dotorminn  thr  trHnr;mittPd  frequency  to  +  2  opa .  The  nn~ 
xuro  of  tdie  meaGurod  doppler  data  in  sufficiently  smooth  that  there  is  no 
difficulty  in  matching  curves  to  this  accuracy. 


Computation  of  Doppler  Frequency  Shift 

The  radio  doppler  frequency  shift  caused  by  a  moving  transmitter  is 


given  by 


Af  ^  _  n  dr 
fo  c  dt 


where  fo  is  the  transmitter  frequency 

and  r  is  the  range  along  the  curvdd  ray 

path.  dr  _  q  dE  .  n  ^  /rt\ 

^  -  p  cos  p  ^  sin  p  ^  (8) 

where  E  is  the  angle  at  th^  earth* s  center, 
between  the  satellite  and  the  receiver. 


\  / 
'Ve-  E  -W 


Applying  Snel 1  * s  lav'  for  spherical  siratification 
n  a  cos  a  -  np  cos  B 
we  have  finally: 

dr  _  no  rdE  ,  t'^nB  dpi 

—  ~  —  a  cos  I—  +  — ^ 
d  t  n  n  Cl  t  p  cit 

Simil ar ly 


p  cos  (B  +  6^)  ^  +  sin  (p  +  6p) 


2'  dt 


where  R  is  the  straight  line  range. 

Combining  (?) ,  (9),  and  (lO)  gives: 

^  ^  -no  j-cos,^  ^  __  g_co3_qo  ^  (p  ^  ^  ).tan  p)]  (]]) 

fo  c  ‘‘cos  a  dt  p  dt  '  2' 

If  both  6  and  6^  are  small,  (11)  becomes: 
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(12) 


2  2 

Af  .  -no  r/n  6  .  ,  N  dR  n  s  p  tan  bn  dpn 

"  T  ^  a)  ^  a 

o 

In  cur  computations  we  actually  use 

^  .  >1  ^  / 

fo  c  d  t 

The  errors  due  to  refraction  (difference  between  (12)  and  (13)  are  at 


most  1  cpG  at  100  me  when  a  -  0. 

^  0 

dR 

Tc  conipuLe  —  of  tne  satellite,  we  proceed  as  follows; 
-  2ap  cos  E 

R  d£/a  r£  _  „-|  _£  d_co^.E  ,  ^  \ 


a  dt 


satellite 


In  the  figure,  0  is  tne  observing  point  at  latitiide  ==  ^ 

the  orbit  inclination, 

Cos  E  -  A  sin  D  t  B  cos  D 

where  A  -  sin  I  sin  lat  t  cos  I  cos  lat  sin  AI^ 

o  0 

and  B  -  cos  lat  cos  AT 
0 

d  cos  E  _  /*  ^  ^  dD  .  .  .  dAL  r  .  ^ 

- T~: —  ~  (A  cos  D  B  sin  D)  tt  cos  lat  — rr  Lsin  D  cos  i  cos  AL 

at  dt  o  dt 

-  cos  D  sin  AT-] 

If  6  Is  the  angular  distance  of  the  satellite  from  perigee  and  r  is  the 
length  of  the  semimajor  axis,  and  if  the  orbit  follows  Kepler's  Laws, 


d  p/a  d6  r (l  -  e  )  e  sin  9 


dt  dt 


j-U  -  e  ;  e  sin  o 
(lie  cos  0)^ 


Furthermore 


d0  ^  2>.' 

~  Y  ^  ^  (1  -  e  ) 
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Where  T  i3  the  satellite  period  and  e  is  the  orbit  eccentricity. 


d  AL  ^ 
dt 


Also  '"jj.  -*  earth  rotation  rate  +  time  rate  of  change  of  the  R.  A.  ‘^f 


ascending  node. 

It  js  more  convenient  to  change  from  0  to  a  variable  cp,  defined  as 


3ln  e  (3  - 


e+  cpsG 


in  cp  =  - - n  and  cos  cp  -  ■7-7 - IT — n 

]  +  e  cos  6  ^  1  +  e  cos  U 


with  this  substitution,  (14)  becomes  finally 

~d^  ”  R  ~  ^  cp  -  (]  -  e^)^(A  cos  D  -  B  sin  1))] 

C15) 

-  ^  cos  ]at  [sin  D  cos  I  cos  AL  -  cos  D  sin  AL  ]| 

Q  0  dt  0  0  j 

In  this  expression  D  is  actually  0  +  co  (argument  of  perigee)  and 
AL  is  Long  Aries  -  R.  A.  of  ascending  node  -  Longitude  of  observer. 

The  change  to  the  variable  cp  gives 

Y  “  (t  -  tp)  =  cp  -  e  sin  cp 
where  tp  is  time  of  perigee. 


Computational  Procedure 

From  (15)  it  will  be  noted  that  the  dominant  term  in  the  one  involving 
A  cos  D~B  sin  D,  "which  is  zero  at  tan  D  =  A/B.  This  corresponds  closely  to 
the  point  of  zero  doppler  shift. 

It  was  assumed  that  the  least  well  known  quantity  in  (15)  is  D.  For 
this  reason,  the  time  of  nearest  approach  is  determined  approximately  from 
the  experimental  doppler  data  (time  of  greatest  rate  of  doppler  change). 

The  value  of  D  is  computed  at  this  time  as  tan  D  =  A/B  where  A  end  B  ere 
determined  from  given  orbital  elements.  This  gives  reference  values  of 
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0«  tp  and  t  “  tp.  Computations  are  performed  at  this  time  and  at  equally 
spaced  time  intervals  on  either  side.  For  each  time  interval  change  cp  is 
deterTTined  by 


A  cp 


m 


m  e  sin  cp  ^ 


m 


ecoscp^-^  -  2(]-e  cos  9)  2(]-ecoscp) 


r  2  ,2 

Le  sin  cp 


]  -  e  cos  cp 


e  cos 


<p]) 


where  m  -  constant 


The  new  value  of  cp  “  9^  +  /Ap  is  then  used  to  compute  the  doppler  at  the 
new  time..  Since  the  right  ascension  of  the  ascending  node  and  the  argu¬ 
ment  of  peiigee  change  by  at  most  a  few  degrees  per  day,  it  is  assumed  that 
these  /alues  are  constant  during  a  20  minute  pass. 


Comparison  with  Observed  Doppler  Data 

In  practice,  it  is  found  that  the  observed  curve  of  doppler  frequency 
versus  time  rarely  matenes  the  computed  curve  within  the  t  1  cps  tolerance 
allowed  by  tne  measured  data.  Differences  of  the  order  of  +  50  cps  are 
usually  found..  In  order  to  bring  the  curves  into  coincidence,  it  is  neces¬ 
sary  to  modify  tne  given  orbital  elements  by  a  small  amount.  In  the  case  of 
Transit  II,  the  Applied  Pnysics  Laboratory  of  the  Johns  Hopkins  University 
re~det ermines  the  orbital  elements  each  day.  It  has  been  found  that  even 
these  elements  require  small  corrections  in  order  to  apply  to  the  portions 
f  the  paSvSes  wnich  are  visible  from  our  station.  Corrections  of  the  nodal 
right  ascen.sinn  of  the  order  of  a  few  tenths  of  a  degree  are  usually  suf¬ 
ficient  to  match  the  doppler  data.  Occasionally  it  is  necessary  to  also 
move  the  argument  of  perigee  by'  a  degree  or  two.  Using  these  relatively 
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small  changes  the  doppler  computations  of  most  Transit  passes  have  been 
brought  into  coincidence  with  the  observed  data  to  within  +  1  cps . 

In  the  case  of  Tiros  I,  the  situation  is  somewhat  different.  Since 
Che  eccentricity  is  much  lower  than  for  Transit,  (.0045  compared  to  ,031) 
reasonable  changes  in  eccentricity  and  the  argument  of  perigee  are  not 
effective  in  changing  the  doppler  curve.  However,  the  assumption  of  con¬ 
stant  transmitter  frequency  during  the  pass,  which  is  required  for  a  valid 
doppler  match,  is  not  met  nearly  as  well  with  Tiros  as  with  Transit. 

This  is  reasonable,  since  great  pains  were  taken  to  design  the  Transit 
system  with  this  in  mind,  whereas  the  Tiros  transmitter  was  not  so  designed. 
Therefore  the  limit  of  accuracy  with  which  one  can  match  doppler  curves  is 
not  83  good  as  the  doppler  measuring  system  will  allow.  On  some  passes, 
good  matches  are  obtained  with  slight  changes  in  nodal  right  ascension. 

On  ether  passes,  it  is  evident  that  the  transmitter  frequency  is  drifting 
by  a  few  cps  daring  the  pass,  such  that  one  cannot  match  the  entire  pass 
with  any  one  set  of  reasonable  orbit  elements.  Since  one  has  no  way  of 
knowing  the  way  in  which  the  transmitter  frequency  is  changing,  one  is  less 
certain  of  the  exact  orbital  elements  which  apply.  On  the  Tiros  passes 
which  are  considered  to  be  valid,  the  match  is  of  the  order  of  +  5  cps.  In 
most  of  these  cases  the  uncertainty  in  the  true  elevation  angle  trajectory 
lb  of  the  Older  of  +  0,':)  mr  or  less.  During  half  of  the  Tiros  passes  (6 
=  ’Ut  of  13)  w'hich  were  otherwise  suitable  for  analysis,  no  reasonable  doppler 
match  could  be  made,  indicating  a  frequency  drift  during  the  pass  of  at 
least  10  cp".  Because  of  the  uncertainty  of  the  true  trajectories,  these 
passes  were  discarded.  This  situation  occurred  on  only  one  of  the  14 
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Transit  passes  during  July,  when  the  transmitter  was  operating  properly. 
During  August  however,  difficulties  with  the  solar  batteries  deteriorated 
the  Transit  frequency  stability  to  the  point  where  only  one  ou;^  of  three 
otherwise  suitable  passes  could  be  well  matched. 

Frc'm  the  above  discussion  one  can  conclude  that  the  satellite  trans¬ 
mitter  frequency  must  be  constant  during  a  15  minute  interval  to  +  1 

g 

part  in  10  in  order  to  allow  computation  of  the  true  trajectory  to  +  0.1  mr. 
Measured  Elevation  Angle  Error 

Figure  9  shows  an  idealized  plot  of  the  apparent  elevation  angle 
versus  time,  measured  with  the  sea  interferometer.  The  shaded  areas 
represent  the  ’’ smeared-out  "  regions  within  which  the  source  lies,  due  to 
/cry  rapid  position  fluctuations,  and  the  center  dots  represent  the  measured 
positions  corresponding  to  passage  through  the  minima.  It  can  be  seen  that 
the  points  show  second- to-second  variations,  but  the  set  follows  a  well 
defined  trend.  The  solid  line  in  the  figure  represents  the  true  elevation 
angle  eirov  a  plotted  versus  time.  The  elevation  angle  error  b  caused  by 
atmospheric  refraction  of  the  radio  wave,  is  simply  the  difference  between 
a  given  point,  and  the  ordinate  of  the  curve  at  the  same  time.  The  dashed 
curve  represents  the  predicted  trajectory,  based  on  corrections  computed 
as  described  in  Section  VI.  Figures  10  through  37  are  plots  of  these 
measured  values  of  b  as  a  function  of  the  apparent  elevation  angle  a^. 

In  thes^'  plots,  the  ordinate  of  scale  is  ten  times  the  abscissa  scale. 

It  can.  be  seen  that  all  of  these  figures  show  a  similar  trend  of  decreasing 
6  with  increasing  a^.  In  some  of  the  figures,  a  distinct  difference  in 
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refraction  is  evident  between  the  two  antenna  heights.  This  is  attri¬ 
butable  to  the  fact  that  on  these  days  a  strong  inversion  layer  was  known 
to  exist  at  a  height  about  even  with  the  upper  antenna.  In  this  situation, 
the  sea  reflected  wave  is  bent  an  additional  amount  in  passing  through 
the  invf^rsion  layer  while  the  direct  wave  to  the  upper  antenna  Is  not. 

The  principal  effect  is  that  the  reflected  wave  arriving  at  the  antenna 
has  its  reflection  point  moved  a  little  closer  to  the  antenna,  resulting 
in  a  greater  path  difference  for  a  given  elevation  angle  of  the  incoming 
wave.  Thus,  the  elevation  angle  where  a  given  minimum  will  occur  is  a 
little  less  than  if  thelayer  were  not  present.  The  direct  wave  to  the 
lower  antenna  on  the  other  hand  will  also  be  bent  by  the  inversion  layer, 
resulting  In  valid  apparent  elevation  angles  in  this  caseo  Computation 
shows  that  for  a  strong  inversion  layer  (AN  -  60  unit  s)  at  the  antenna 
height,  this  effect  will  result  in  elevation  angles  which  are  about  2.1  mr 

too  low  at  a  -  0,  0.6  mr  at  100  mr  and  0.3  mr  at  200  mr .  Since  these  are 
0 

the  order  of  magnitude  of  the  observed  differences,  as  well  as  the  observed 
AN’s,  it  is  evident  that  one  should  use  the  values  measured  at  the  lower 
antenna  in  these  cases.  For  comparison  purposes  of  course,  one  must  compute 
the  6’r  with  the  invei sion  layer  above  the  antenna,  as  will  be  discussed 
later . 

Effects  of  Sea  Swell 

An^'lher  feature  of  the  6  plots  (Figures  10  through  37)  is  the  presence 
of  a  periodic  fluctuation,  particularly  at  higher  elevation  angles.  It 
has  been  found  that  the  time  periods  involved  are  of  the  order  of  15  seconds 
which  corresponds  to  that  of  sea  swells.  The  most  pronounced  effect  occurs 
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cn  passes  urhere  the  satellite  was  to  the  southwest  during  the  interfero¬ 
meter  portion  of  the  pass,  although  a  less  distinct  effect  is  also  observed 
n  westerly  azimuths,  with  least  effect  on  northwesterly  azimuths.  The 
swells  were  observed  to  be  aligned  along  a  SW-NE  orientation,  which  would 
make  them  aligned  along  the  look  direction  in  the  most  pronounced  cases, 
and  normal  to  the  look  direction  in  the  least  pronounced  cases.  The  distance 
between  crests  for  a  15  second  period  is  of  the  order  of  1000  feet. 

Computation  of  the  size  of  the  first  Fresnel  zone  gives  dimensions  of 
the  order  of  0.1  mile  widtn  by  one  mile  length-  It  is  clear  that  when  the 
"  ng  dimension  of  the  zone  is  aligned  with  the  swell  orientation,  the 
average  level  of  the  water  surface  will  rise  along  its  entire  length  as 
the  swell  crest  passes  by,  since  the  Fresnel  zone  is  not  wide  enough  to 
encompass  both  the  crest  and  trough  simultaneously. 

When  the  look  azimuth  is  normal  to  the  swells  on  the  other  hand, 
several  crests  and  troughs  vrill  be  illuminated  simultaneously,  resulting 
in  a  periodic  residual  effect  which  is  much  smaller  than  in  the  former  case. 

The  effect  of  sea  swells  will  be  greatest  at  high  elevation  angles. 
Since  a  --  n\/2ri,  it  follows  that  Aa  alAh/h),  At  200  mr  for  example,  a 
2  foot  ii^e  in  effective  sea  level  will  cause  a  lowering  of  the  minima 
pattern  by  0  5  mr  ,  and  at  50  m-^  ,  the  effect  would  be  about  0.1  mr.  The 
bserved  periodic  fluctuations  are  in  agreement  with  these  figures. 

The  shape  of  sea  swells  is  characterized  by  long  flat  troughs  with 
relatively  short  peaked  crests.  Therefore,  the  average  sea  level  will 
correspond  more  nearly  to  the  level  of  the  troughs  tnan  with  the  crests. 

It  would  appear,  therefore,  that  the  computea  elevation  angles  of  the 
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minima  correspond  most  closely  with  the  lower  water  levels  (higher  angle 
values).  Thus  the  true  measured  6’s  are  just  below  the  peaks  of  the  periodic 
fluctuations  found . 


In  the  top  section  of  Table  I  the  representative  measured  6’ 
tabulated  for  each  pass,  at  of  150  mr,  100  mr ,  50  mr,  20  mr , 
10  mr.  These  data,  together  with  the  remainder  of  the  table  will 
cussed  in  Section  VII „ 


are 

and 

be  dis- 


-  37  - 


VI.  COMPUTATION  OF  6  VS  a  FROM  ROUTINE  REFRACTIVE  INDEX  DATA 

0 

The  calculation  of  refractive  effects  over  a  spherical  earth  with  a 
spherically  stratified  auaosphere  involves  knowing  the  vertical  distribu¬ 
tion  of  refracijve  index-  Since  this  distribution  varies  considerable 
both  geographically  and  in  time,  it  is  necessary  that  the  computation  method 
be  able  to  handle  any  given  distribution.  The  method  developed  by  Smyth 
Research  Associates  handles  this  problem  in  a  particularly  simple  way- 
The  atmosphere  is  divided  into  a  series  of  layers,  each  of  which  has  a 
linear  change  of  refractive  index  with  height  (constant  gradient)- 

In  such  a  layer  it  can  be  sho^m  that  the  refractive  bending  A  y 

exactl’y.  ^ 

Ay  2iNg-  N^)  X  ]0  V(tan  +  tan 

where  N  is  the  atmospheric  refract ivity 

and  j3  15  the  angle  which  the  ray  makes  with  the  horizontal  at 
a  given  level. 

The  pubscripls  B  and  T  refer  to  the  bottom  and  top  of  the  layer, 
respectively- 

Similarly  it  can  oe  shown  that  the  range  error,  caused  by  retardation 
-f  the  ray  relative  to  vacuum  velocity,  is  very  nearly; 

-  (Mg  +  N^)  (H^  -  Hg)  X  Pg  +  sin  p^) 

where  the  'H’ s  are  the  heights  of  the  top  and  bottom  of  the  layer- 
The  ierivati'n  of  these  expressions  is  straightforward,  though  somewhat 
tedi'^Uo.  The  reader  is  referred  to  the’ fol  1  owing  paper  for  details: 

"Simple  Method  for  C  mputing  Tropospheric  and  Ionospheric  Refractive 
Effects  on  Radio  Wave^",  by  S.  Veisbrod  and  L-  J.  Anderson,  Proceedings 
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47,  pp.  ] 770-3777,  October,  3959. 


Trr~^Dosphere 

T  compute  tropospheric  effects,  one  uses  radiosonde  data.  These 
data  are  in  the  form  of  temperature,  relative  humidity,  and  pressure  at 
what  are  termed  ’'significant  levels".  These  levels  have  been  selected 
during  analysis  of  the  original  traces,  such  that  between  any  two  adjacent 
levels,  the  temperature  and/or  humidity  varies  very  nearly  linearly. 

These  levels  then  define  the  boundaries  between  layers  of  constant  refrac¬ 
tive  index  gradient,  as  required  by  the  expressions  given  above.  The 
radiosonde  data  at  each  level  are  converted  into  refractivity  N,  which  is 
(n  -  1)  X  10^  (where  n  is  the  actual  refractive  index)  by  the  following 
expression . 

N  =  (77.6/T)  (P  t  4810  e/T) 

w^here  T  is  the  temperature  in  degrees  absolute 

P  is  the  atmospheric  pressure  in  millibars 
and  e  is  the  vapor  pressure  in  millibars. 

T  complete  the  conversion,  (  ne  also  needs  to  determine  the  height  at 
Which  each  significant  level  occurs.  This  information  is  normally  supplied 
by  the  meteorologist  who  works  up  the  original  radiosonde  trace,  or  else 
It  can  be  readily  '^btained  from  available  tables. 

We  now  have  N  values  and  corresponding  heights  at  each  significant 
level.  Next,  tan  P  is  determined  at  each  level,  using  Snell’s  law: 

tan^P  [(l  +  -)(l  -(n  -  N,  )  X  10  cos  a  ]  ^ 

nan  o 

One  determines  tan  P  at,  each  significant  level  for  each  of  interest. 
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Using  these  values  and  the  N  values  at  corresponding  heights  enables  one  to 
c'mput'-e  Ay  for  each  layer  whose  boundaries  are  adjacent  significant  levels. 
Finally,  one  sums  up  thcAy^s  for  all  layers  to  obtain  total  y^  for  the 
troposphere . 

I  :'nosphei  e 

We  must  now  compute  similar  values  of  y^,  the  ionospheric  ray  bending. 
The  refractive  index  data  we  need  for  this  purpose  are  in  the  form  of 
icnogram  traces,  which  can  be  converted  to  electron  density  versus 
height.  The  ionospheric  refractivity  corresponding  to  a  given  electron 
density  is  also  a  function  of  the  radar  frequency  (unlike  the  tropospheric 
N^s  w^hich  are  independent  of  frequency,) 

The  expression  for  ionospheric  refractivity  is: 

N  ■■  -4.03  N  X  ]0"Vf^ 

1  e 

wnere  is  electron  density  in  el ectrons/meter^ 

and  f  is  the  radai  frequency  in  megacycles/second. 

The  ionosphere  is  divided  into  layers  of  constant  N  gradient  as  be¬ 
fore  ana  the  ray  bending  is  again  computed  as  t.bove . 

HeVing  determined  y^  at  each  ,  we  now  add  the  tropospheric  and 
in'opheric  contributions  to  obtain: 

Y  Tt  ^  Ti 

Since  we  are  interested  in  the  elevation  angle  correction  6  rather 
than  the  ray  bending  itself,  we  must  correct  for  the  fact  that  the  trans¬ 
mitter  (• r  radar  target)  is  not  at  infinite  distance,  where  b  and  y  are 
equal.  We  determine  b  as  follows; 
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-■  (y  tan  p  +  Y^/2  ^  ^)/(y  +  tan  P  -  tan  a^) 

In  this  case,  the  P’s  refer  to  the  height  of  the  transmitter  or  radar 
taiget.  They  are  computed  from  Snell’s  Lav  given  in  the  tropospheric 
section . 

Ran^e  Error 

Tne  radar  range  error  in  the  troposphere  is  computed  according  to  the 
expression; 

AR  (ft.)  -  (Ng  +  N^)  (h^  -  h^)/lC00  (sin  p^  +  sin  p^) 

where  the  h’s  are  expressed  in  thousands  of  feet. 

Again  the  ZXR’s  are  computed  for  each  layer  and  added  to  give  the  to¬ 
tal  AR  fox  the  troposphere. 

I  nospheric  retardation  is  computed  using  the  N^’s  obtained  from  the 
i^n'gram  trace.  Again  the  AR‘s  are  added  to  give  total  ionospheric  re¬ 
tardation.  The  total  retardation  for  the  entire  atmosphere  is,  of  course, 
the  sum  of  tropospheric  and  ionospheric  contributions, 

D  pT)]er  Error 

The  d^ppler  frequency  shift  which  one  measures  after  a  radio  v/ave 
has  traversed  the  atmosphere  differs  from  what  it  would  be  for  a  transmitter 
-r  radar  target  moving  in  exactly  the  saiTie  way,  in  the  absence  of  an  atmos¬ 
phere.  In  presen^je  of  refractive  bending,  the  measured  doppler  corresponds 
t‘  a  wave  leaving  the  transmitter  or'  target  along  the  ray  path,  whereas  in 
absence  of  refraction,  the  wave  would  leave  along  the  straight  line  path  to 
the  receiver.  This  slight  difference  in  angle  is  equal  to  y  “  Therefore, 
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the  doppJer  error  is: 

Af  =  (y  -  6)  AF 

where  /SF  is  the  measured  doppler  shift.  Since  y  -  b  rarely  ex¬ 
ceeds  one  milliradian  at  satellite  or  missile  heights,  it  follows  that  the 
dnpp; er  error  rarely  exceeds  0 :  f  the  measured  doppler. 

Measured  dopplers  are  of  the  order  of  a  few  kilicycles,  thin,  the 
doppler  errors  caused  by  refraction  are  of  the  order  of  a  few  cycles  at 
most . 
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COMPARISON  OF  MEASURED  AND  PREDICTED  ELEVATION  ANGLE  ERROR 


In  table  I,  the  measured  values  of  elevation  angle  error  6  are 
tabulated  at  apparent  elevation  angles  of  150,  100,  50,  20  and  10 
mil  1 iradians .  In  the  next  lower  part  of  the  Table,  the  tropospheric  6^s 
computed  f'^om  radio  sonde  data  taken  at  Pt .  Arguello  on  the  coast  50  miles 
southwest  of  the  station,  arc  tabulated.  Two  radiosonde  ascents  were 
made  each  day,  and  the  resulting  6*s  were  computed  for  the  above  angles 
for  each  ascent.  These  were  then  interpolated  in  time  for  comparison  with 
measured  passes. 

The  next  lower  part  of  the  table  shows  the  remaining  6,  which  is  the 
ionospheric  component.  The  succeeding  lines  of  the  Table  show  the  predicted 
ionospheric  component,  based  upon  four  objective  methods. 

The  first  of  these  involves  using 'the  NBS  Series  D  Charts,  “Basic 
Radio  Propageilion  Predictions."  These  charts  allow  one  to  determine  the 
n  nthly  average  vertica''  MI^F  of  the  ionosphere  corresponding  to  the  point 
of  F  layer  penetration  Dy  the  "^adio  wave  on  each  satellite  pass.  Using 
i hese  data  and  the  method  described  in  SRA  Report  91 ,  one  can  estimate 
the  expected  ionospheric  component  of  5  for  each  frequency. 

i'he  second  set  of  predictions  utilizes  ionogram  data  measured  at 
White  Sands,  New  Mexico,  and  Boulder,  Colorado,  on  the  days  the  passes 
occurred.  Such  data  are  taken  at  15  minute  intervals,  and  the  ionograms 
measured  closest  to  the  same  local  time  as  tnit  corresponding  to  ray 
penetration  of  the  ionosphere  on  each  pass,  were  used  to  compute  indivi¬ 
dual  Ionospheric  refraction  ermponents.  The  assumption  which  is  made  is 
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TABLE  I 


COMPARISON  OF  MEASURED  AND  PREDICTED  EI.EVATION  ANGLE  ERROR 


Transit  Par.ses 


No. 

2 

II 

12 

15 

16 

17 

25 

Yve(^ 

.  (me) 

216 

162  216 

162  216 

216 

162  216 

162  216 

162 

216 

o(  (inr) 

n  ’ 

Measured  h  in  milliradians 

150 

5.9  3.8 

4.1 

7.1  4.1 

5.1  .*.6 

5.7 

4.3 

1  oc 

3.4 

4.7 

7.6  5.2 

4.6 

8,0  5.0 

6.0  5.3 

6.7 

5.2 

50 

5.6 

6 . 5 

10.5  8.0 

6.8 

1 — 1 

o 

8.1  7.2 

9.7 

8.1 

20 

8.1 

11.0  10.2 

16.4  13.1 

10.3 

1 — 1 

00 

CO 

1 — 1 

12.4  11.3 

14.7 

13.4 

10 

10.1 

13.2 

14.0 

16.8 

Tropospheric  ^  Computed 

from  Radiosonde  Data 

1  50 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2 

.0 

IOC 

2.9 

2.9 

2.8 

2.8 

2.8 

2.8 

2 

Q 

•  / 

5C 

5.C 

5.0 

4.7 

4.9 

5.0 

5.0 

5 

X 

or' 

< 

9.7 

8.1 

8 . 6 

8.8 

8.8 

8.8 

ip 

^ "  .8 

12.9 

12.6 

13.2 

Ionospheric  Component 

(Mea 

surer]  minus  Trcposoheric) 

1  r,' 

3.9  1.8 

2  J 

5.1  2.T 

3.1  2.6 

3.7 

2,3 

KX 

0.5 

1  .8 

4.8  2.4 

.8 

5.2  2.2 

3.2  2.5 

3.8 

2.3 

50 

0.6 

T  5 

5.8  3.3 

1  .9 

5.3  2.6 

3.,1  2.2 

4.7 

3.1 

2^ 

-0.6 

4.2  1.4 

8.3  5.0 

1  .7 

IC.O  4.6 

3.6  2.5 

5.9 

4.1 

10 

-1  .'7 

0.3 

1  .4 

3.6 

Predicted 

Ionospheric 

Compenents  (Average  10-3  50  inr 

) 

m3 

Series  D 

1  .0 

1.8  1.0 

1.9  1.0 

1  .0 

1.9  1.0 

1.7  0.9 

1.5 

0.8 

I  •' n 

gram 

0,8  0.5 

0.8  0.5 

0.5 

2.1  1.2 

1,9  1.1 

1.1 

0.6 

Mear 

c<.  (iiff 

’ .8  1  .0 

5.3  3.0 

6.9  3.9 

1.6  0.9 

3.4 

1  .9 

Me.i. 

0.5 

3.2  1  .8 

3.6  2.0 

2.1 

4.6  2.6 

3.6  2.0 

3.8 

2.2 
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TABTE  I  (Cont'  l) 


Vac  3 

Nc'. 

2£' 

28 (e fit) 

29 

30 

31  (ri 

se) 

31 (eel 

Freq 

.  (me) 

162 

216 

162  216 

162  216 

162  216 

162 

216 

216 

T.r) 

Mea.sured 

^  in  Mill iradians 

< 

O 

4.6 

2.9 

5.2 

5.2  3.5 

5-3  4.7 

100 

6.2 

4.2 

6.4 

6.0  4.4 

6..]  5.5 

5.4 

50 

9.2 

7.2 

12.1  9.1 

8.7  6.8 

7.6  6.7 

7.3 

6.2 

8.3 

20 

13.7 

1]  .4 

14.8 

12.8  ]]  .0 

J1.3  10.0 

11.8 

10.7 

13.7 

iO 

14.8  13.9 

13.9 

Tropo 

spheric  ^  Computed  from 

Radiosonde 

Data 

15C 

2 

,  1,  ' 

2.0 

1  .9 

1  .8 

1  . 

8 

1  .8 

1  Qn 

2 

.9 

2.9 

2.8 

2.7 

2. 

6 

2.6 

50 

e 

.r 

5.0 

4.8 

4.7 

4. 

5 

4.5 

20 

P 

.p 

8.8 

8.2 

8.0 

n 

1  . 

5 

7.5 

]C 

11  .9 

]0, 

,0 

.  I'n  'Gpliei  iC  C-'r"poncnt  .'j'leasured  minus  Tiopospneri 

■c) 

’  50 

2 « 6 

0.9 

3.2 

3.3  1.6 

3.5  2.9 

5.3 

.1  .3 

1.5 

3.2  1.6 

3.4  2.8 

2.8 

5'J 

4.2 

2.2 

7.1  4.1 

3.9  2.0 

2.9  2.0 

2.8 

1  .7 

3.8 

^ 't 

4.^ 

2.6 

6.0 

4.6  2.8 

3.3  2.0 

4.3 

3.2 

6.2 

r 

2.9  2.0 

3.9 

Predicted  I -nr^spheric  Components  (A/erage 

10-150 

mr) 

Serlf^r 

■  .6 

0.9 

2.7  1.5 

1.6  0.9 

1.8  1.0 

1  .5 

0.8 

1  .5 

i  n  ' 

rra^. 

’  .  ’’ 

<1.6 

1.7  "  .0 

1.1  0.6 

1.1  0.6 

1  .1 

0.6 

Mea. 

.iiff 

.4  .  3 

2.4 

6.'i  3.9 

4.0  2.3 

2.1  1.2 

2.5 

1.4 

Me  a:: 

2.2 

'  .2 

5.7  3.2 

3 .  1.7 

4-0  2.3 

2.8 

1  .6 

2.8 

_  /  _ 


TABLE  I  (Cont'd) 


Pass  No, 

32 

33 

36(ri3e) 

36(set) 

Freq.  (me) 

162 

216 

362  216 

162  216 

324  362 

216 

o(^(mr) 

150 

5.6 

4.3 

Measured 

4.2  3.4 

oC  in  Milliradians 

0 

5.7  4.5  4.9 

3.1 

TOO 

6.4 

4.7 

5.2  4.6 

6.7  5.7 

4.9  6.0 

4.1 

50 

8.5 

6.5 

7.2  6.8 

9.3  8.3 

7.4  7.9 

6.1 

20 

13.0 

3  1  .2 

33  .2  30.5 

32.8  11.8 

9.2 

10 

15.5  34.6 

35.5  35.0 

12.6 

Tropospheric  6  Computed  from  Radiosonde  Data 

150 

1 

.9 

3  .9 

2.0 

2 

.0 

100 

2 

.8 

2.8 

2.9 

2, 

.9 

50 

4 

.9 

4.9 

5.0 

5, 

.0 

20 

8 

.4 

8.4 

8.7 

8 

.7 

10 

32.2 

12.7 
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TABLE  I  (Cont'd) 


Tiros  Passes  (]08  me) 


Pass  No. 

] 

3 

5 

40 

41 

43 

47R 

47S 

57 

59R 

59s 

61 

(mr) 

Measured 

b  in  Milliradians 

150 

5.3 

4.1 

3.3 

4.9 

9.3 

TOO 

6.6 

6.7 

6.4 

6 . 1 

4.4 

7.4 

10.5 

8.7 

6.0 

50 

8 .6 

8.9 

8.5 

9.5 

7.8 

11 .0 

10.3 

13.6 

12.6 

9.3 

11  .2 

c  2 

20 

11.2 

13.3 

11.7 

12.1 

15.8 

15.3 

16.6 

16.1 

14.3 

17.4 

I4.I 

10 

12.9 

13.4 

22.3 

21 .8 

19.5 

17.8 

Tropospheric 

b  Computed 

from  Radiosonde  Data 

150 

2.0 

2.0 

2.0 

2.1 

2.1 

2.1 

2.0 

2.0 

2.1 

2.1 

2.1 

2.1 

100 

2.9 

2.9 

2.9 

3.0 

3.0 

3.0 

2.9 

2.9 

3.0 

3.0 

3.0 

3.0 

50 

4.9 

5.0 

5.1 

5.1 

5.1 

5.1 

5.0 

5.0 

5.1 

5.1 

5.1 

5.1 

20 

8.6 

8.7 

8.8 

8.7 

8.7 

8.7 

8.7 

8.7 

9.1 

9.1 

9.1 

8.6 

n  r\ 

J  0 

12.2 

11.5 

11.8 

11  .8 

11.8 

11.3 

lonovjpheric  C 

omponent  (Measured 

minus 

Tropospheric) 

^50 

3.3 

2.0 

]  .2 

2.9 

7.3 

100 

3.7 

3.8 

3.5 

3.1 

1.4 

4.5 

7.6 

5.7 

3.0 

50 

3.7 

3.9 

3.4 

4.4 

3.7 

5.9 

5.3 

8 .6 

7.5 

4.2 

6.1 

4.1 

20 

2.6 

4 .6 

2.9 

3.4 

7.1 

6.6 

7.9 

7.0 

5.2 

8.3 

^  .5 

10 

0.7 

1 .9 

10.5 

10.0 

7.7 

6.5 

Predicted 

Ionospheric  Component 

(Average  10- 

■1  50  mr) 

iIB3  Sene..  D 

3 .6 

3.6 

3.6 

1  .8 

2.1 

3.2 

4.7 

7.2 

4.3 

4.9 

6.8 

2.1 

I  nogram 

0.9 

]  .] 

2.5 

6.6 

8.0 

2.9 

3.6 

1.4 
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that  the  ionospheric  structure  is  only  ?j  function  of  local  time  (solar 
elevation  angle)  and  that  other  random  changes  over  a  two  hour  period  are 
neg] igibl e . 

The  geomagnetic  latitude  of  White  Sands  is  the  same  as  Pt.  Buchon, 
and  that  of  Boulder  is  10*^  further  north.  Thus  on  nearly  westerly  look 
angles,  the  White  Sands  data  were  used,  and  on  n'^rthwester'ty  angles,  Boulder 
data  were  used.  There  are  no  nearby  ionosonde  staxions  whose  geomganetic 
1  at  nude  iu  about  10°  south  of  Pt,  Buchon;  south  and  southwesterly  passes 
are  therefore  apt  to  be  underestimated,  since  electron  density  is  usually 
higher  at  lower  latitudes. 

The  detailed  ionospheric  computations  were  made  in  accordance  with 
pr-^cedures  outlined  in  the  previous  section. 

The  thiid  orediction  method  requires  measured  data  at  two  fre¬ 
quencies.  Fi  om  the  measured  differences  in  at  these  frequencies,  and 
The  ascomption  that  the  ionospheric  b  is  inversely  proportional  to  the 
frequency  squared,  one  can  determine  the  ionospheric  6  at  each  frequency. 

The  fourth  method  requires  measured  6's  at  two  frequencies,  and  is  in 

reality  a  medoure  of  the  self-consistency  of  all  passes  and  of  the  ^ 

f 

aossumption 

In  Table  11  the  average  residual  error  for  all  passes  where  the  re¬ 
quired  data  were  available,  is  tabulated  at  several  elevation  angles,  for 
eacn  uf  the  fouT  methods  described  above.  The  numbers  in  parentheses  are 
averages  of  less  than  6  passes,  and  are  therefore  not  reliable.  It  can 
be  seen  that  all  of  the  prediction  methods  give  residual  errors  which  usu¬ 
ally  increase  vitn  decreasing  frequency  and  with  decreasing  The 
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TABLE  II 

AVERAGE  RESIDUAL  ERROR  IN  5(mr.) 


f  (me 

)/oC^I,mr)  10  20 

50 

100 

3  50 

Radiosonde 

and 

Scries  D 

108 

3.5  3.7 

3  .3 

0.6 

0.7 

162 

(3  .3)  2.9 

2.2 

3 .8 

3  .6 

216 

3.7  2.2 

3  .3 

3  .2 

3  .2 

Radiosonde  and 

,  lonogram 

108 

(4.2)  2.4 

2.4 

3  ,5 

3.4 

■'62 

1 — 1 

2.5 

1  .9 

3  .8 

2]  6 

3.5  2.3 

3  .5 

3  .2 

3  .2 

'62 

Measured 

(3.3)  ".6 

oC  Differences 
0 

0.8  ].0  1.2 

2l6 

1.5  1.3 

0.8 

3  .0 

3  .3 

M.ea.sured 

High 

Ang3  e 

6’ s 

TC8 

3.’’  1.7 

1  .3 

0.7 

0.0 

"62 

(3  .0)  "  .7 

0.8 

0.4 

0.3 

216 

i — I 

0.5 

0.3 

0.2 

two-frequency  measured  differences,  yield  residual  errors  which  are 
appreciably  leas  than  the  Series  D  method,  and  which  are  reasonably  the 
same  at  both  frequencies.  Finally,  the  self-consistency  test,  labelled 
"Measured  High  Angle  b’s",  shows  residual  errors  which  are  much  smaller 
at  high  angles,  but  about  the  same  at  low  angles. 

From  the  data  indicated  in  the  lower  sections  of  Table  II,  one  can 

conclude  that  the  ^  assiamption  is  well  substantiated  by  the  measurements, 
f^ 

Furthermore,  the  self-consistency  of  the  two-frequency  (Transit)  passes 
is  good  (+  0.3  mr)  at  high  angles  and  the  results  increase  as  the  eleva¬ 
tion  angle  decreases.  This  indicates  that  the  tropospheric  component 
predictions  are  underestimated,  since  ionospheric  6's  are  almost  constant 
over  the  measuring  interval  . 

It  is  further  apparent  that  the  individual  ionogram  based  predictions 
are  no  better  than  those  based  on  monthly  averages  of  the  NBS  Seried  D 
circular,  and  are  if  anything  somewhat  worse.  The  only  conclusion  one  can 
iraw  frcm  this  comparioon  is  that  the  assumption  negligible  random  ionos¬ 
pheric  structure  changes  over  a  two  hour  interval  is  not  as  valid  as  the 
assumption  that  our  passes  occurred  during  conditions  which  were  close  to 
the  monthly  average.  Another  factor  which  undoubtedly  contributed  to  the 
underestimation  of  the  ionogram  6’s  is  the  previously  mentioned  lack  of 
lata  for  southerly  passes. 

It  appears  that  the  only  way  in  which  to  improve  the  ionospheric 
predictions  would  he  to  have  ionogram  data  taken  near  the  actual  region 
Where  the  ray  paths  penetrate  the  F  layer.  Since  this  takes  place  in 
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a  ring  between  750  and  ]]00  miles  from  the  station  and  at  azimuths  from 
160  to  325,  thivS  would  be  difficult  to  accomplish  by  vertical  incidence 
soundings.  However,  the  u.ie  of  the  oblique  incidence  technique  from  che 
station  at  Stanford  would  provide  ionospheric  information  from  the  actual 
penetration  region.  Since  tne  times  and  look  dii actions  of  satellite  passes 
can  be  predicted  with  adequate  accuracy  tor  two  v/eeks  or  more,  it  would 
be  possible  to  schedule  times  and  azimuths  for  the  oblique  incidence 
soundings  well  in  advance,  to  insure  obtaining  the  required  ionospheric 
data,  without  the  necessity  for  a  blanket  coverage  during  the  expected 
days  of  pass  observation.  It  is  suggested  that  this  approach  be  investi¬ 
gated,  since  it  represents  about  the  best  one  can  hope  to  do  in  obtaining 
ionospheric  information  for  refraction  predictions. 

As  mentioned  earlier,  the  tropospheric  component  was  computed  from 
ladiosonae  data  taken  twi.e  a  da;y  at  a  point  50  miles  away  from  the  re- 
oeiving  jite  ,  Inaomucn  as  the  ray  paths  in  the  crucial  lowest  30,000 
feet  of  tno  atmosphere  occupy  about  200  miles  in  length,  it  would  seem  that 
the  50  mi'^es  separation  would  not  be  significants  Since  both  locations 
a’^e  on  the  coastline,  the  measured  profile  would  be  as  representative  of  the 
ray  path  whetner  tnc  ua^a  were  taken  at  Ft.  Aiguello  or  at  Pt ,  Buchon. 

One  adv/antage  may  wel^  ^Couli  f  i  urn  a  sounding  taken  near  the  actual  time 
'f  tne  measu'.ed  pass.  Whether  or  not  tne  hour-to-hour  changes  measured 
at  the  Lite  would  be  reflected  in  similar'  changes  throughout  the  important 
tropospheric  part  of  the  ray  patru  ,  remains  to  be  seen.  Radiosonde  equip¬ 
ment  is  invj tailed  at  the  field  site  for  obtaining  these  data,  but  during 


tne  period  of  the  meaGurements  reported  herein,  was  not  in  operating 
condition.  Future  pas^^es  which  are  measured  will  be  accompanied  by  on-site 
radiosonae  data  in  order  to  test  this  point. 

The  resulto  ...nown  in  Table  1 1  indicate  that  the  radiosonde-derived 
tropospheric  oorreflions  are  generally  too  small  at  low  elevation  angles. 

This  is  typical  of  strong  ’"efrac  .ion  conditions. 

It  was  realized  that  the  radiosonde  sensing  elements  (particularly 
the  humidity  strip)  have  appreciable  lag  coefficients,  and  therefore 
tend  to  indicate  inversion  layers  which  are  thicker  and  more  diffuse  than 
Lhey  actually  are.  This  benavior  would  resul t  in  6*s  which  are  too  small 
at  low  elevation  angles,  as  is  evident  in  Tables  I  and  II.  In  reducing 
tne(  radiosonde  data,  therefore,  compensation  was  made  for  the  lag  in  the 
humidity  strip.  It  was  assumed  tiiat  above  the  base  of  the  inversion  l^er 
the  vapor  piOw^sare  decreased  as  rapidly  as  the  temperature  increased,  and 
reached  its  equllibiium  value  at  tne  top  of  the  inversion  layer.  The 
v«por  presrure  value  used  was  the  equilibrium  value  at  the  top  of  the 
inversion  layer.  The  vapor  pressu'^e  value  used  was  the  equilibrium 
value  whicn  was  indicated  at  a  neight  500  to  1,000  feet  above  the  top  of 
the  inversion  layer.  Il  is  felt  that  most  of  the  effects  of  strip  lag 
have  been  eliminated  oy  this  procedure.  The  lag  of  the  temperature  element 
at  tne  dustomary  rise  rate  of  1000  ft/minute  is  about  5  seconds..  Its  effect 
IS  to  lower  the  in.^ersicn  layer  by  100  feet  at  most,  which  has  no  appreciable 
effect  on  the  resulting  6  values. 

Figure  jB  1C  a  plot  of  tne  computed  6'g  versus  No,  and  indicates  the 
variability  caused  by  the  shape  of  the  individual  N  profiles.  It  can  be 
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seen  that  at  200  mi,  the  computed  6^0  average  about  .02  mr  above  the 
standard  atmosphere  with  a  variability  of  +  .01  mir .  At  100  mr,  the 
corresponding  figures  are  +  .03  and  +  .02  mr  respectively.  Even  at  an 
elevation  angle  as  low  as  50  mr,  the  actual  profile  shapes  yielded  average 
5^s  only  0,1  mr  above  the  standard  atmosphere  with  a  variability  of  +  .05 
mr .  kt  angles  of  50  mr  and  higher,  therefore,  one  could  conclude  that  the 
effects  of  the  strong  refraction  conditions  encountered  are  quite  small, 
and  can  be  well  accounted  for  by  standard  atmosphere  corrections  with  a 
small  systematic  increase;  or  by  corrections  based  upon  No  alone. 

At  10  mr  and  20  mr,  the  layer  effects  become  progressively  more  pro¬ 
nounced.,  The  detailed  6  computations  give  values  generally  from  0.5  to 
1  mr  higher  than  the  standard  atmosphere  at  20  rr.r  and  1  to  3  mr  higher 
at  10  mr.  The  pi ofil e-to-profil e  variability  is  +  0.2  mr  at  20  mr  and 
+  0 .6  mi  at  10  mr . 

It  should  be  realized  however,  that  these  compai isons  are  between 
standard  atmospnere  profiles  and  the  moderate-to-strong  refractive  pro¬ 
files  measured  by  the  Pt.  Arguello  radiosonde  (corrected  for  lag  as  des¬ 
cribed  above) .  The  Table  II  data  indicate  the  corresponding  additional 
differences  between  ob.served  o's  and  those  shown  in  Figure  38. 

Thus  at  10  and  20  mr  ,  while  the  computed  6’s  are  about  2.0  and  0.7 
mr  higher  than  standard  respectively,  the  observed  6*s  are  an  additional 
1 .4  mr  higher.  At  50  mr  we  must  add  about  0.8  mr  to  the  0.1  mr  computed 
excess  over  standard.  At  higher  angles  the  average  residual  tropospheric 
corrections  are  about  10  times  the  profile-minus-standard  differences,. 
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In  summary  then,  one  can  conclude  that  comparisons  such  as  that  shown 
in  Figure  38  tend  to  give  one  a  false  sense  of  security  about  the  accuracy 
with  which  corrections  can  be  made  for  tropospheric  refraction*  The  actual 
uncertainty  at  angles  above  50  mr  is  of  the  order  of  10  times  as  much  as 
one  would  expect  on  the  basis  of  N  profile  computations  alone.  At  lower 
angles  the  ratio  decreases  to  about  two  to  one,  but  the  absolute  value  of 
the  uncertainty  about  doubles  in  going  from  an  of  50  mr  to  10  mr. 

These  comparison  v;ere  made,  of  course,  under  strong  to  moderately 
strong  refraction  conditions  and  as  such,  represent  about  the  worst  un¬ 
certainties  one  can  expect  to  enccunter . 
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VIII.  SHORT  TERM  FLUCTUATIONS 


Elevation  Angle  Error 

Examination  of  Figures  10  to  37  indicates  that  the  second-to-second 
fluctuations  in  6  at  108  me  are  of  the  order  of  +  1  mr  at  elevation  angles 
from  0  -  20  mr  and  +  0.5  nir  at  50  mr .  Corresponding  values  at  216  me  are 
0.4  mr  and  0.2  mr.  The  presence  of  sea-swell  fluctuations  at  higher  angles 
makes  it  difficult  to  estimate  the  fluctuations  at  greater  than  50  mr, 

but  it  is  certain  that  they  will  decrease  rapidly  with  increasing  At 

^  '3  greater  than  150  mr ,  the  fluctuations  are  probably  smaller  than  the 
0."'  mr  resolution  of  the  measuring  system. 

One  can,  however,  draw  certain  conclusions  about  very  rapid  fluctura- 
tions  in  6  by  examining  the  max-min  ratio  of  the  interference  patterns.  If 
one  had  a  point  source  of  radiation  and  a  smooth  reflecting  surface  with 
unity  divergence  factor,  the  depth  of  individual  minima  would  be  limited 
only  by  the  receiver  noise  level.  At  low  elevation  angles,  however,  the 
divergence  factor  is  less  than  unity  and  increases  as  increases.  The 
effect  of  wave  height  on  the  amplitude  of  the  reflected  component  on  the 
'ther  hand  is  small  at  low  angles  and  increases  with  Figure  39  shows 

the  greatest  max-rcin  ratio  achievable  as  a  function  of  for  various 
ratios  of  wave  height  to  radio  wavelength.  At  108  me  (\  =  9.1  Tt) ,  a  wave 
height  of  4.5  ft.  v;ould  limit  the  ratio  to  l6  db  for  the  lowest  minimum, 
but  the  ratio  increases  to  a  maximum  of  32  db  at  -  37  mr.  Measured 
max-min  ratios  of  as  much  as  30  db  are  frequently  observed  although  ratios 
as  little  as  18  db  are  also  seen. 

If,  instead  of  being  a  point  source,  the  apparent  elevation  angle  of 
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the  satellite  fluctuated  very  rapidly  (with  millisecond  periods)  the  source 
would  appear  to  have  a  finite  angular  width  In  this  case  the  max-min 

ratio  would  be  limited  by  A^,  and  one  should  be  able  to  draw  cohclusions 

about  the  maximum  value  of  from  the  observed  ratios. 

One  can  show  that  under  ideal  reflection  conditions  the  max-min  ratio 

would  be.  ^  ^  ('jThA)^(cos  ^1^-2  sin  /\X)^ 

0  o 

Where  h  is  the  height  of  the  antenna  above  the  reflecting  surface,  \  is 
the  wavelength,  and  is  the  wavelength  spread  caused  by  the  radio  band¬ 
width.  In  our  system,  we  can  use  h  ==  860  ft.  and  a  bandwidth  of  500  cps. 

This  gives  A\  5  x  10  '  ft.  at  100  me  and  1.2x10  £1  at  216  me.  At 

-  200  mr  ,  tp(  of  1  mr  would  thus  correspond  to  an  R  of  5,2  db  at  2l6  me 
and  11  .2  db  at  106  me,  and  a  A^  of  0.1  mr  would  correspond  to  an  R  of  25.0 
db  at  21 6  me,  and  29-8  db  at  108  me.  Thus  the  minimum  observed  R  of  about 
18  db  coirespnds  to  a  Ax  rf  the  order  of  0.2  mr  at  2l6  me  and  0.5  mr  at  108 
me.  The  maximum  R  rf  about  30  db  corresponds  to  a  "'f  0.05  mr  at  2l6 
me  and  0.1  mr  at  108  me.  The  fact  that  R*s  of  as  much  as  30  db  have  been 
observed  indicates  that  A^X  cannot  be  greatei  than  the  minimum  values  cited 
above.  If  it  were  greater  than  this,  f ' r  a  given  pass  no  R  would  be  as 
gre.at  as  30  db.  The  fact  that  some  R’s  less  than  this  are  observed  can 
be  attributed  to  fluctuations  is  sea  roughness  and  other  causes  not  related 
tne  atmosphere,.  It  is  therefore  concluded  that  the  common  occurrence 
of  3C  db  max-min  ratios  indicates  a  maximum  millisecond  fluctuation  ampli¬ 
tude  in  6  order  of  0.^  mr  and  0,05  mr  at  108  me  and  216  me  respectively 
even  at  elevation  angles  as  low  as  30  mr. 


The  doppler  data,  which  were  oi'igina]]y  taken  for  radio  fix  purposes, 
are  also  a  very  sensitive  means  of  determining  range  eiror  fluctuation 
rates.  If  one  imagines  a  wave  penetrating  a  section  of  the  ionosphere 
with  a  given  integrated  electron  cnunt  (between  the  transmitter  and  the 
receiver),  the  retardation  of  the  wave  will  be  proportional  to  the  inte¬ 
grated  count..  If  now  the  ray  path  sweeps  through  an  ionospheric  blob 
having  a  higher  electron  density,  the  retardation  will  increase,  causing 
a  decrease  in  doppler  frequency  which  would  be  proportional  to  the  gradient 
in  retardation  normal  to  the  ray  path.  For  example,  if  the  range  error 
were  3  km  during  a  given  one  second  interval ,  and  if  it  increased  by  1% 
during  the  next  one  second  interval,  the  apparent  range  change  rate  would 

be  1%  of  3  km  =  30  meters  per  second.  This  would  correspond  to  a  doppler 
30  -7 

ohift.  of  — -  10  ,  or  10  cp^.  at  100  me.  Since  the  doppler  caused  by 

3x10^ 

motion  f  the  satellite  is  a  smooth  function  of  time,  one  can  easily 
distinguish  the  second-to-second  departures  caused  by  electron  density 
gradients  ^rnm  the  doppler  caused  by  satellite  motion. 

Figures  40  and  41  f'lustrate  observed  doppler  fluctuations  on  two 
passes  of  Tiros  and  Transit  respectively.  Several  other  such  passes  have 
been  plotted,  and  the  results  are  similar  to  those  shown.  It  is  seen  that 
during  mOot  of  the  pass  w^hich  is  above  the  horizon,  the  average  fluctuation 
IS  of  the  order  cf  +  2  cps  at  108  me  and  +  *1  cps  at  216  me. 

The  data  plotted  in  Figures  40  and  41  are  actually  the  differences 
between  the  received  frequency  averaged  over  a  given  one  second  interval 
and  that  averaged  over  a  succeeding  one  second  interval  starting  a  second 


beyond  the  first  internal ,  Thus  one  has  averaged  readings  spaced  two  sec¬ 
onds  apart.  Each  pass  is  plotted  in  two  parts  in  order  to  spread  out  the 
time  scale.  The  pass  starts  in  each  case  at  the  left  ordinate,  near  the 
zero ^ difference  line,  and  proceeds  across  the  width  of  the  figure,.  The 
remainder  of  the  pass  starts  again  at  the  left  ordinate  (at  a  value  of 
-10  to  -20)  and  continues  across  the  figure.  The  average  fluctuation 
rates  snown,  inviicate  rates  of  change  of  range  error  of  the  order  of  3 
meters  per  second  at  108  me  and  0.8  meters  per  second  at  216  me. 

In  general,  the  fluctuations  are  reasonably  random,  although  one 
occasionally  finds  several  consecutive  intervals  during  which  the  doppler 
rate  is  consistently  below  average  or  above  average.  This  would  indicate 
that  the  range  erroi  probably  fluctuates  about  the  mean  value  with  as 
much  as  +30  meter  and  9  meter  deviations  at  108  me  and  216  me  respectively 
over  pe:'-iods  of  10  seconds. 

These  are  occasional  intervals  (Figure  41 )  where  the  doppler  fluc- 
tuationo  exceed  5  cps  per  second,  indicating  deviations  of  the  order  of 
30  meters  in  a  two  ,:econd  interval  .  In  one  case  (just  before  the  geo¬ 
metric  norizen  point,  Pass  16,  Figure  41 )  a  cumulative  range  deviation 
)f  12  meters  was  reachea  before  it  swung  back. 

As  pointed  o-jr,  above,  these  range  eiror  deviations  are  visualized 
as  being  caused  by  the  ray  patn  sweeping  through  blobs  with  higher  or 
'’eiiser  electr  n  density  than  their  surroundings.  It  appears  from  the 
observed  doppler  data  nowever,  that,  when  averaged  over  one  second  inter- 
va'ls,  the  ionosphere  13  remarkably  smooth.  This  means  either  that  the 
electron  density  gradients  take  place  over  distances  of  many  tens  of  miles 
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Guch  that  the  gradientG  themselveG  are  low,  or  that  there  is  a  very  fine 
grain  structure  (less  than  a  mile  in  extent)  which  irS  averaged  out  during 
the  one  second  counting  interval  . 

There  wil"^,  of  coux’se,  always  be  some  averaging  because  of  the 
finite  length  of  the  Fresnel  zone,  which  at  108  me  is  of  the  order  of  6 
miles  by  one  mile  in  dimensions  at  low  elevation  angles,  3.7  miles  by 
one  mile  at  300  mr,  and  one  mile  by  one  mile  at  the  zenith.  These 
dimensions  assume  a  satellite  height  of  500  miles  and  the  ionospheric 
"layer’’  height  at  300  km.  Since  the  satellite  velocity  is  roughly  4-5 
miles  per  second,  one  could  obtain  a  resolution  of  a  mile  by  recording 
doppler  data  at  0.2  second  intervals  instead  of  one  second  intervals. 

At  21 6  me,  the  Fresnel  zone  dimensions  are  rexluced  to  0.7  of  the  above 
va'^ues.  Stfl  higher  frequencies  would  have  smaller  Fresnel  zone  sizes, 
but  wou^d  also  be  less  sensitive  to  electron  density  gradients.  It 
would  appear  that  21 6  me  would  allow  one  to  obtain  optimum  overall  re¬ 
solution  of  range  error  fluctuations,  of  the  order  of  0.7  mi'^es.  Doppler 
resolution  rf  0.2  cps  in  a  time  interval  of  0,2  second  would  allow  one  to 
detect  range  error  rate  change.^  from  mile-to-mile  of  the  order  of  one 
meter  per  second . 
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IX.  PRECURSORS 


One  uf  the  results  of  the  investigations  described  in  this  report 
was  quite  unexpected.  This  was  the  detection  of  signals  many  minutes 
before  the  satellite  rose  above  thehorizon,  or  after  it  had  set.  The 

strength  of  these  "precursors"  is  occasionally  within  a  few  db  of  the 

free  space  levels.  The  more  distant  signals  usually  rise  above  the 
detection  level  for  a  few  seconds  and  then  fade  out  for  a  minute  or  more 
and  recur  again  at  closer  range.  The  signals  do  not  show  any  rapid  fad¬ 
ing  structure  (within  the  .01  second  response  cf  the  recorder) .  Occa¬ 
sional  ly  there  is  a  slower  fading  period  of  a  second  or  so  with  fades 
of  the  order  of  a  few  db,  suggestive  of  a  multipath  mechanism. 

Out  of  a  total  of  4.7  passes  during  which  these  signals  would  have 
been  detected,  they  actually  appeared  on  35.  In  Figure  42,  the  times  of 

lise  or  set  are  plotted  versus  azimuth.  The  purpose  of  this  plot  was  to 

determine  whether  the  pattern  would  be  as  one  would  expect  for  an 
ionospheric  mechanism.  Characteristically,  the  MIJF  f  the  ionosphere 
decreases  north  of  the  station  latitude  (35^N)  and  increases  to  the  south. 
It  is  also  at  a  minimum  at  midnight.  If  the  ionosphere  were  responsible 
therefore,  one  would  expect  the  "no-show"  passes  to  occur  toward  the  north 
■I  at  night.  The  data  plotted  ^n  Figure  42  suggests  such  a  pattern,  but. 
there  are  some  inconsistencies  where  a  precursor  was  detected  at  a 
particular  time  and  azimuth  during  one  day  and  not  on  the  next.  This 
w"uld  be  expected  if,  as  the  measured  b  data  indicate,  the  horizontal 
pattern  of  MUF  varies  considerably  from  day  to  day. 
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Tn  Figure  43,  the  level  of  the  precursor  peaks  is  plotted  versus  the 
distance  of  the  satellite  telov  the  radio  horizon  of  the  receiver.  This 
distance  was  computed  along  the  line  connecting  the  satellite  and  re¬ 
ceiver  locations  at  the  time  of  precursor  occurrence.  The  ordinate  is 
expreesed  in  db  below  the) free  space  level  at  a  distance  corresponding  to 
the  first  two  or  three  maxima  in  the  interferometer  pattern.  The  free  space 
■’evel  to  this  point  cou*^!  be  established  since  at  the  maximum  the  signal 
should  have  been  about  5  db  above  free  space  because  the  direct  and 
reflected  waves  add  In  phase  at  these  points. 

It  can  be  seen  that  the  most  distant  signals  were  received  when  the 
satellite  was  nearly  2,000  mi'^es  beyond  the  radio  horizon,  and  that  the 
signal  levels  were  10  and  18  db  below  the  free  space  level  at  this  dis¬ 
tance.  Several  other  signals  were  detected  at  1200  mile  ranges,  including 
^ne  case  just  2  db  below  free  space. 

In  examining  these  data,  no  well  defined  frequency  dependence  is 
evident.  Whi'^e  most  of  the  long  range  signals  were  at  the  higher  fre¬ 
quencies,  some  of  them  show  162  me  stronger  than  21 6  me  and  others  are 
reversed.  If  the  ionosphere  were  responsible  as  the  propagation  mechanism, 
one  would  expect  the  108  me  signals  to  be  much  more  predominant  than  those 
at  higher  frequencies. 

Correlation  of  the  dates  of  precursor  occurrences  with  radiosonde 
data  suggests  that  tropospheric  ducts  are  probably  responsible  for  the  ob¬ 
served  signals.  On  the  days  when  the  strongest  inversion  layers  were 
present,  the  precursors  were  strongest  and  most  persistent.  During 
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days  of  relatively  weak  inversion  layers,  the  precursors  tended  to  be 
weak  or  absent. 

Since  most  of  the  observations  were  made  in  the  presence  of  inversion 
'ayers,  it  Is  desirable  to  extend  them  to  a  wider  variety  of  conditions 
in  order  to  observe,  for  example,  whether  precursors  ever  occur  in  the 
absence  of  such  layers.  Once  the  precursor  mechanism  has  been  positively 
identified,  it  will  be  possible  to  estimate  the  early  detection  possibi¬ 
lities  of  this  phenomenon  in  various  areas  of  operational  importance. 
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X .  CONCLUSIONS 


The  work  which  has  been  described  in  this  report  has  led  to  several 
conclusions  regarding  the  accuracy  with  which  one  can  correct  for  atmos¬ 
pheric  refractive  effects  on  precision  location  of  objects  above  the 
ionosphere . 

1.  The  average  total  elevation  angle  error  measured  during  this 
work  is  shown  in  the  following  table: 


TABLE  III 
Average  6*s 


freqy^ 

3-0 

20 

50 

]00 

150 

108 

]8.0 

^4.4 

]0.0 

7.0 

5.4 

162 

15.3 

13.3 

9.0 

6.4 

5.4 

2]6 

]3.8 

T]  .4 

7.2 

4.9 

4.0 

?.  The  correction  'f  elevation  angle  errors,  using  nearby  radio¬ 
sonde  data  and  NBS  monthly  ionospheric  predictions  as  well  as 
daily  ionograin  data,  leaves  residual  errors  averaging  +  1-3 
milliradians  at  angles  above  3°.  Corresponding  values  at  1^ 
are  2.5  mr . 

3.  Much  of  this  residual  error  is  caused  by  space  and  time  variations 
in  the  ionosphere,  which  are  not  derivable  from  the  monthly 
average  ionospheric  predictions  or  from  daily  ionogram  data  taken 
1000  miles  away.  Measured  differences  in  elevation  angle  at  two 
frequencies  allow  one  to  evaluate  ionospheric  corrections  which 
reduce  the  residual  errors  appreciably  from  the  values  given 
above . 
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4.  Further  analysis  snows  that  the  remaining  part  of  the  residual 
errors  at  low  angles  (l  to  3  degrees)  is  attributable  to  under¬ 
estimation  of  tropospneric  refraction  in  spite  of  compensating 
for  a  lag  of  the  radiosonde  sensing  elements. 

3.  The  self  cnnsistenty  of  14  passes  where  two-frequency  data  were 
available  indicates  residual  errors  of  0.3  milliradians  at  angles 
above  5^. 

6.  Second-tO“3econd  fluctuations  in  apparent  elevation  angle  at 
108  me  average  about  1  mr  at  angles  below  1°  and  about  0o5  mr 

at  3^.  The  216  me  fluctuations  are  about  half  the  108  me  values. 

7.  An  estimate  of  very  rapid  elevation  angle  fluctuations  (periods 
v':>f  tlie  order  of  milliseconds)  indicates  a  maximum  amplitude  of 
0.'^  mr  at  108  me  and  0.05  mr  at  216  me  This  estimate  is  based 
on  the  observed  max-min  ratios  of  the  interferometer  patterns. 

8.  The  3econd-tO“second  fluctuations  in  range  error  can  be  estimated 
from  measured  doppler  fluctuations.  These  are  usually  small; 

'f  tne  order  of  +  3  meters  at  108  me  and  0.8  meters  at  216  me. 
Occasional  deviations  of  ten  times  these  values  are  found o  De¬ 
viations  greater  tnan  12  meters  at  216  me  are  rare. 

0,  Precursor  signals  nave  been  received  on  10%  of  the  passes  where 
they  would  have  been  detected.  These  signals  were  received  when 
the  satellite  was  as  much  as  1800  miles  beyond  the  horizon.  No 
well  defined  frequency  dependence  is  observed  on  the  occurrence 
or  intensluy  of  the  signals.  Passes  where  no  precursors  were  found, 
were  generally  tnose  requiring  a  northwesterly  look  angle,  or 
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those  occurring  early  in  morning.  The  propagation  mechanism 
producing  precursor  signals  is  ^ery  efficient;  signals  of  nearly 
free  space  intensity  were  observed  out  to  distances  of  1200  miles 
beyond  the  radio  horizon.  Correlation  with  radiosonde  data  indi¬ 
cates  that  the  mechanism  responsible  for  propagating  precursor 
signals  is  the  topospheric  duct  formed  by  the  elevateded  inver¬ 
sion  layer  which  was  present  during  the  passes  when  strong  signal 
were  observed . 
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XI.  RECOMMENDATIONS 


As  a  result  of  the  information  obtained  to  date  the  following 
recommendationwS  are  made: 

^ .  Obtain  oblique  incidence  ionuspjiere  sounding  data  which  are 

simultaneous  with,  and  from  the  same  area  as  F  layer  penetration 
of  the  ray  paths  for  each  measured  pass.  The  station  at  Stanford 
would  be  best  suited  for  the  Pt,  Buchon  meavsurements . 

2.  Obtain  on-site  radiosonde  data  siniul taneous  with  each  pass,  for 
computation  of  the  tropospheric  refraction  components. 

3.  Further  emphasis  should  be  placed  on  using  the  doppler  fluctu¬ 
ation  approach  for  obtaining  finer  grained  information  on  ionospheric 
inhomegeneities . 

4.  Further  investigation  f  precursor  signals  should  be  performed, 
particularly  during  weak  duct  and  standard  tropospheric  conditions 
to  verify  whether  or  not  tropospheric  ducting  is  the  only  cause 

of  such  signals. 

5.  Further  measurements  of  elevation  angle  errors  should  be  made 
during  weak  duct  and  standard  conditions,  so  the  conclusions 
presented  above  may  be  extended  to  a  wider  are  a  of  applicability. 
Winter  and  springtime  measurements  at  Pt.  Buchon  wculd  provide 
such  information.  Coverage  of  a  wide  variety  of  conditions  would 
enable  extropol ation  of  position  accuracy  limitations  to  any  area 
of  operational  interest. 

6.  Simultaneous  refraction  measurements  should  be  made  using  radio 
stars  and  radio  satellites.  This  would  indicate  the  utility  of 
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radio  star  monitoring  as  a  means  of  determining  day-to-day  refrac¬ 
tive  corrections  on  an  operational  basis.  Such  measurements 
should  be  made  over  an  annual  cycle  to  provide  the  full  diurnal 
range  of  ionospheric  conditions  at  the  time  of  setting  of  the 
star . 
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